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SPECTRAL ASYMPTOTICS FOR V- VARIABLE SIERPINSKI 

GASKETS 

U. FREIBERG, B.M. HAMBLY, AND JOHN E. HUTCHINSON 


Abstract. The family of V -variable fractals provides a means of interpolating 
between two families of random fractals previously considered in the literature; 
scale irregular fractals (V = 1) and random recursive fractals (V = oo). We 
consider a class of V -variable affine nested fractals based on the Sierpinski 
gasket with a general class of measures. We calculate the spectral exponent 
for a general measure and find the spectral dimension for these fractals. We 
show that the spectral properties and on-diagonal heat kernel estimates for 
V -variable fractals are closer to those of scale irregular fractals, in that it is 
the fluctuations in scale that determine their behaviour but that there are also 
effects of the spatial variability. 


1. Introduction 

The field of analysis on fractals has been primarily concerned with the construc¬ 
tion and analysis of Laplace operators on self-similar sets. This has yielded a well 
developed theory for post critically finite (or p.c.f.) self-similar sets, a class of 
finitely ramified fractals [31]. One motivation for the development of such a theory, 
aside from its intrinsic mathematical interest, has come from the study of transport 
in disordered media. However, in this setting the fractals arise naturally in mod¬ 
els from statistical physics at or near a phase transition and are therefore random 
objects without exact self-similarity but with some statistical self-similarity. 

In order to develop the mathematical tools to tackle analysis on such random 
fractals one approach has been to work with simple models based on self-similar sets 
but exhibiting randomness. The first case to be treated was that of scale irregular 
fractals [17], [2], [24] and [11], which have spatial homogeneity but randomness 
in their scaling. A more natural setting is provided by random recursive fractals, 
initially constructed by [38], [12], [16], where the fractal can be decomposed into a 
random number of independent scaled copies. The study of some analytic properties 
of classes of random recursive Sierpinski gasket can be found in [18] , [20] and [22] . 

Recently there has been work tackling random sets arising from critical phenom¬ 
ena directly, with a particular focus on the percolation model. Substantial progress 
has been made in the study of random walk on critical percolation clusters in the 
high dimensional case, see [3] and [36]. A bridge between these two approaches 
can be found in work on the continuum random tree [9], [10] or on critical perco¬ 
lation clusters on hierarchical lattices [23], both of which have random self-similar 
decompositions and hence have descriptions as random recursive fractals. 

In this paper we consider V-variable fractals recently introduced in [6, 7]. This 
class of random fractals is defined via a family of iterated function systems and 
a positive integer parameter V. It interpolates between the class of homogeneous 
(scale irregular) random fractals, corresponding to V = 1, and the class of random 
recursive fractals, corresponding to V = oo. As for the random recursive fractals 
we can regard these V-variable fractals as determined by a probability measure on 
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the set of labelled trees. In this case the measure is not a product measure, but is 
defined in a natural (if not completely obvious) manner which allows for at most 
V distinct subtrees rooted at each level. 

Our aim in this paper is to investigate the analytic properties of the class of 
R-variable Sierpinski gaskets and to compare their behaviour to the scale irregular 
and random recursive cases. We show their Hausdorff dimension in the resistance 
metric is the zero of a certain pressure function and their spectral dimension, the 
exponent for the growth of the eigenvalue counting function, is the zero of another 
pressure function. The connection between these two dimensions is established. 
We develop and extend standard methodology to examine more detailed properties 
of the eigenvalue counting function and the on-diagonal heat kernel. These results 
show that the R-variable fractals are closer to the scale irregular case, in that 
their fine properties are generally determined by fluctuations in scale rather than 
fluctuations which occur spatially across the fractal. 

Model problems. We consider two model problems. Recall from [26] the descrip¬ 
tion of a self-similar set as an iterated function system (or IFS) at each node of a 
tree generated by the address space. 

Homogeneous and Random Recursive Fractals. For the first model problem we con¬ 
sider the two IFSs generating the Sierpinski gasket fractal SG(2) and the fractal 
SG(3) defined in [17]. The scale factors for SG(2) are mass m 2 = 3, length £2 = 2 
and time S 2 = 5. For SG(3) we have mass m 3 = 6 , length £3 = 3 and time 
S 3 = 90/7. The conductance scale factors can be computed directly, or from the 
Einstein relation p = s/m, giving p 2 = 5/3, /O 3 = 15/7. Let (M,S,L) be a triple of 
random variables taking each of the values (m.j, Si,£/) where * = 2, 3 with probabil¬ 
ities p, 1 — p respectively. 

Then, for the V = 1 (homogeneous) case, we construct a random fractal using a 
sequence taking its values in {2,3} and applying the corresponding IFS to all sets 
at a given level of construction. A realization of the first few stages can be seen 
in Figure 2. Then a simple scaling analysis shows that the Hausdorff dimension 
is given by df = E log M /E log L where E denotes the expectation with respect to 
the probability measure generating the sequence. For the spectral dimension with 
respect to the natural “flat measure” one can extend the idea from [14] and [35] 
in the case of a single IFS fractal and apply a scaling argument to the Dirichlet 
form together with a Dirichlet-Neumann bracketing argument, see [19]. This gives 
the spectral dimension d s = 2E log M/E log S. For the V = 00 (random recursive) 
case, each IFS is chosen independently for each node at each level. In this case 
we have d s = 2 d r ^/(d r j + 1) where d r f is the Hausdorff dimension in the resistance 
metric, that is d r j is such that E (M(S/M)~ d f) = E M 1+dr f S~ d f = 1. The argument 
again uses scaling properties of the Dirichlet form and a Dirichlet-Neumann brack¬ 
eting argument, see [19, 20]. An alternative approach to computing the spectral 
dimension for random V = l,oo fractals is via heat kernel estimates, see [ 2 ] and 
[17, 18, 19, 20], 

The second model problem is drawn from the class of affine nested fractals con¬ 
sidered in [13]. This model interpolates between the slit triangle (which is not itself 
an affine nested fractal) and SG(3). Consider 7 triangles in the configuration shown 
in Figure 1 and take £ as the side length of the three triangles at the corners of 
the original triangle. The side lengths of the other triangles are given as 1 — 2t for 
the three triangles on the centre of each side and 31? — 1 for the downward pointing 
central triangle, where 1/3 < £ < 1/2. As 1! 1/2 we have the slit triangle and at 

£ = 1/3 we have SG(3). We construct a homogeneous random or random recursive 
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Figure 1. A member of the family of Sierpinski gaskets interpo¬ 
lating SG(3) and the slit triangle, where 1/3 < £ < 1/2. 


fractal by taking a suitable distribution for £ on [1/3,1/2) and either using a se¬ 
quence, applying the same IFS at each node in the construction tree for the F = 1 
case, or independently for each node in the V = oo case. 

We note that even scale irregular (F = 1) affine nested gaskets of this type 
have not been treated before and as a consequence of our results we will be able 
to calculate the Hausdorff and spectral dimension for the random homogeneous 
version (V = 1). By the triangle-star transform, if we assume that the resis¬ 
tance of each piece is proportional to its length, then the resistance scale factor is 
(2£ + l)/(£ + 2) in that if we take resistances on the three different types of trian¬ 
gle to be (£ + 2)/(2£ + 1)(£, 1 — 2£, 3£ — 1) then this is electrically equivalent to the 
triangle with unit resistance on each edge. 

In Section 2 we recall from [8] the Hausdorff dimension result for F-variable 
fractals, and we derive the spectral dimension from our calculations in Sections 4 
and 5. 

V- Variable Fractals. To understand the F-variable versions of our model problems, 
first consider the V = 1 (spatially homogeneous, scale irregular) case of a F-variable 
labelled tree in a manner parallel to the approach taken in the general setting. 
See Figure 2. For F = 1 all subtrees rooted at each fixed level are the same, 
as are the corresponding subfractals at each fixed level, hence the terminology 
“homogeneous”. The subtrees at one level are typically not the same as the subtrees 
at another level, hence the terminology “scale irregular”. 




Figure 2. The level 3 approximation to a 1-variable tree, and 
the prefractal approximation to the associated 1-variable, or scale 
irregular, fractal. Here the family of IFSs is F = {F(2),F(3)} 
with members generating the sets SG( 2) and SG( 3) respectively. 

For a general F-variable tree and for the corresponding F-variable fractal, there 
are at most F distinct subtrees up to isomorphism rooted at each fixed level, and 
correspondingly at most F distinct subfractals up to rescaling at each fixed level 
of refinement. See Figure 3 for a level 2 approximation to a F-variable tree with 
F > 2. In Section 2 we discuss this in some detail and see that there is a natural 
probability distribution on the class of F-variable fractals for each fixed F. 
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The construction of H-variable trees and hence ^-variable fractals will require 
an assignment of a type chosen from {1,..., V}, as well as an IFS, to each node 
of the tree. Nodes with the same type and at the same level will have identical 
subtrees rooted at those nodes. The subfractals corresponding to those nodes will 
be identical up to scaling. See Figure 4. We choose the IFSs according to a 
probability measure and will write Py for the probability measure on the space of 
trees or ^-variable fractals and Ey for expectation with respect to Py. 

Let n( 1) be a random variable denoting the first level after level 0 at which all 
nodes are assigned the same type. Since the number of types is finite and we will 
assume a uniform upper bound on the branching number, Eyn{ 1) < oo. Note that 
n(l) = 1 \iV = 1, and clearly Eyn{ 1) increases with V. 

We write i = i\.. .ik for a node in the tree and denote its height or length by 
|i| = k. The root node is denoted by 0 and |0| = 0. The Hausdorff dimension df 
of the Invariable gasket formed from SG(2) and SG(3) is given Py almost surely 
by the zero of a pressure function in that ( Py almost surely) it is the unique df 
such that Ey logX!|i|=n(i)(^i ‘ ' ^i„(i)) d/ = 0, where £ ik is the length scale value 

1/2 or 1/3 according to which of SG(2) or SG(3) is chosen. See Theorem 2.17, also 
Theorem 2.18. 

Results. For further detail see the Overview at the beginning of the following Sec¬ 
tions 2-5. 

Our main results first establish an expression for the spectral exponent over a 
general class of measures and determine the spectral dimension for these fractals. 
We then provide finer results of two types. We consider the eigenvalue counting 
function and the on-diagonal heat kernel and obtain upper and lower bounds on 
these quantities which hold for all ^-variable trees. By placing a probability mea¬ 
sure on the trees we obtain almost sure results capturing more explicitly their fluc¬ 
tuations. In the model problems the expectation is either over a discrete measure 
on {2,3} or over a suitable distribution on [1/3,1/2]. 

We show in Theorem 4.14 that the spectral exponent can also be expressed as 
the zero of a pressure function. In Theorems 4.16 and 4.18 we see that the spectral 
dimension, the maximum value of the spectral exponent over all measures p defined 
using a product of weights, satisfies the equation d s /2 = d}/(d} + 1) where d} is 
the Hausdorff dimension in the resistance metric. This dimension in turn is the 
zero of another pressure function, see Theorem 3.11. 

We establish upper and lower estimates for the eigenvalue counting function 
and on-diagonal heat kernel for a general class of measures. We show that the 
observed fluctuations arise from two different effects. The first is due to global 
scaling fluctuations as observed for scale irregular nested Sierpinski gaskets [2]. 
The second effect, which arises in the V-variable setting for V > 1 or V = 1 when 
the contraction factors are not all the same, gives additional, though much smaller, 
fluctuations due to the spatial variability of these fractals. 

We first establish from Lemma 4.6 the non-probabilistic result that if Af( A) 
denotes the number of eigenvalues less than A (for the Dirichlet or Neumann Lapla- 
cian), then there is a time scale factor T k: a mass scale factor M k and a correction 
factor A k , such that there are constants c\ , C 2 with 

ciM fe < J\T(A k T k ) and A f(T k ) < c 2 M k , Vfc. 

As in the scale irregular gaskets of [2], this result is true for all realizations. By 
construction the scale factors M k , T k grow exponentially in k but we will be able 
to show that Py almost surely we have A k < ck^ , and even in certain cases A k < 
C(logfc)^, for some constant /?. The spectral exponent for any measure [x defined 
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by a set of weights associated with a given IFS is 

d B (lJ) ,. logA/"(A) 

—— := lim —-—, 

2 K-foo log A 

and we give a formula for this quantity as the zero of a suitable pressure function. 
In the case where the weights are ‘flat’ in the resistance metric we can show that 
there is a function </>(A) = exp(-^/log A log log log A) such that Py-almost surely 

(1) ciA d ' i/2 0(A)- C2 < Af(A) < c 3 A ds / 2 ^(A) C4 , 

for large A, where d s = 2 d r f/(d r f + 1) and d r f is the Hausdorff dimension in the 
resistance metric. 

To compare our results with previous work we note that in the V = 1 case for 
nested Sierpinski gaskets it is shown in [2] that the Weyl limit for the normalized 
counting function does not exist in general and we have for all realizations that 

ci Mk < Af(Tk) < c 2 M k . 


This leads to the same size scale fluctuations as for the V -variable case given in (1). 
For the random recursive case of [20] , the averaging leads to a Weyl limit in that 


lim 

A—>-oo 


N (A) 

\ds/2 


exists Poo a.s., 


where d s = 2 d r j / (dj + 1) and d r j is the Hausdorff dimension in the resistance metric. 

We will also be able to remark on the on-diagonal heat kernel. We note that the 
measures we work with in this setting do not have the volume doubling property 
and hence it is harder work to produce good heat kernel estimates. In the setting 
considered here we can extend the arguments of [2] and [5] to get fluctuation results 
for the heat kernel. In Theorems 5.5 and 5.8 we show that the on-diagonal heat 
kernel estimate is determined by the local environment. In the case where the 
measure is the ‘flat’ measure in the resistance metric we can describe the small 
time global fluctuations in that for almost every point x in the fractal, 

cit~ da,2 (j){l/t)~ C2 < p t (x,x) < c 3 t~ dB/2 (/){l/t) Ci , 0 < t < c 5 , Py a.s., 


for suitable deterministic constants Ci, 02,03,04, and for all t < C5, a random con¬ 
stant depending on the point x. These are of the same order as the V = 1 case 
obtained in [2] and much larger than those in the random recursive case, [22]. 

In the case of general measures we will see that Py-almost surely, /r-almost every 
x in the fractal does not have the same spectral exponent as the counting function 
(except when we choose the flat measure) and thus there will be a multifractal 
structure to the local heat kernel estimates in the same way as observed in [5], [21]. 


We restrict ourselves to affine nested fractals based on the Sierpinski gasket in 
where d > 2. The problem of the existence of a limiting Dirichlet form is not solved 
more generally, even for the case of homogeneous random fractals. If this problem 
were solved, then the techniques used here would enable more general results to be 
obtained concerning V -variable p.c.f. fractals. 


The structure of the paper is as follows. We give the construction of V -variable 
affine nested Sierpinski gaskets in Section 2. We show that by using the structure 
of V -variability there is a natural decomposition of the fractals at ‘necks’; a level 
at which all subtrees are the same. This idea was first used by Scealy in [39]. In 
Section 3 we construct the Dirichlet form, compute the resistance dimension, and 
determine other properties which will facilitate analysis on these sets. In Section 4 
we treat the spectral asymptotics. The heat kernel is dealt with in Section 5. 
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2. Geometry of F-Variable Fractals 

2.1. Overview. Random F-variable fractals are generated from a possibly un¬ 
countable family F of IFSs. Each individual IFS F £ F generates an affine nested 
fractal. We also impose various probability distributions on F. 

For motivation, consider the two model problems in the Introduction. Namely, 
F = { F 2 , F:s} is the pair of IFSs generating SG( 2) and SG( 3), or F is the family 
of affine nested fractals Ft generating the prefractal in Figure 1 for l £ [1/3,1/2]. 

A F-variable tree corresponding to F is a tree with an IFS from F associated to 
each node, a type from the set {1,..., F} associated to each node, and such that if 
two nodes at the same level have the same type, then the corresponding (labelled) 
subtrees rooted at those two nodes are isomorphic. This last requirement is achieved 
by using a sequence of environments , one at each level, to construct a F-variable 
tree. Each F-variable tree generates a F-variable fractal set in the natural way. 
The case F = 1 corresponds to homogeneous fractals and F —> oo corresponds to 
random recursive fractals. 

If all nodes at some level have the same type, the level is called a neck. Neck levels 
are given by a sequence of independent geometric random variables. In Lemma 2.15 
we record some useful results for such random variables. In Section 2.7 we recall 
the Hausdorff dimension result from [8] but in the framework of necks as used in 
this paper, and then give a refinement by using the law of the iterated logarithm. 
This provides motivation for some of the spectral results. 

2.2. Families of Affine Nested Fractals. Let F be a possibly uncountable class 
of IFSs F. each generating a compact fractal K F , and each defined via a set of 
similitudes {iff}i^s F acting on with contraction factors {if } ie ,gF and S F = 
(1,..., N f }. If it is clear from the context we write K , ifi, N and S for K F , iff, 
N F and S F respectively, and similarly for other notation. 

We will have 

3 < iV inf := inf{A F : F £ F} < sup {N F : F £ F} =: N sup < oo, 

^ 0 < := inf{f F : F £ F}. 

The first follows from our later constructions, see (9). The second and third are for 
technical reasons arising in the study of the heat kernel and spectral asymptotics. 
See also the comments after Definition 2.14, from which it is clear that weaker 
conditions will suffice to construct F-variable fractals and establish their Hausdorff 
dimension. 

Let T f denote the set of fixed points of the {iff}i^s F - Then x £ T F is an 
essential fixed point if there exists y £ and i ^ j such that iff (x) = iff ( y ). Let 
Vo denote the set of essential fixed points. 

We always assume that Fo does not depend on F. 

Assume the uniform open set condition for the {iff}- That is, there is a non¬ 
empty, bounded open set O , independent of F, such that {iff ( 0)} ie gF are disjoint 
and U eS F ipf (O) C O. 

Let VC-in = ° ° and let 

(3) V F = |J (Fo), V F =\JV F . 

il,---,i n &S F n>0 

Then K F = cl(Vf ), the closure of V F . 
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For G S F , we call (Vo) an n-cell and ipf t (K F ) an n-complex. 

For x, y £ R d (ir ^ y), set H xy = {z € : |z— x\ = \z —y\} and let U xy : 

be the reflection transformation with respect to H xy . 

When computing the spectral dimensions we further assume each K F is an affine 
nested fractal. That is, the open set condition holds, |Vo| > 2, and: 

(1) K f is connected; 

(2) (Nesting) If (n,--- ,i n ) and (ji, • • • ,j n ) are distinct n-tuples of elements 
from S F , then 

^l-i n (K F ) n^... Jn (K F ) = C-in(^o) (Vo); 

(3) (Symmetry) For x, y € Vq (x ^ y), JJ xy maps n-cells to n-cells, and it maps 
any n-cell which contains elements in both sides of H xy to itself for each 
n > 0. 

We also make the technical assumption that \ip[ (Vo) D iff (Vo)| < 1 for all 1 < 
i^j< N f . 

2.3. Trees and Recursive Fractals. Fix a family F of IFSs as before. For our 
initial purposes it is sufficient only that the IFSs consist of uniformly contractive 
maps on R d . 

Each realisation of a random fractal is built by means of an IFS construction 
tree , or tree for short, defined as follows. 

Definition 2.1. (See Figure 3) An (IFS construction) tree T corresponding to F 
is a tree with the following properties: 

(1) there is a single, level 0, root node 0; 

(2) the branching number N z at each node i has 2 < N 1 < oo (AT* > 3 later); 

(3) the edges with initial node i are numbered (“left to right”) by 1,..., TV*; 
where * = i\ ... ik in the usual manner and |i| := k > 1 is the level of i , or 
i = 0 in which case |z| := 0 is the level; 

(4) there is an IFS F z £ F associated with each node i, TV* = |F*| (the 
cardinality of F z ), and the /cth edge with initial node i is associated with 
the /cth function in the IFS F 1 . 

The unique compact set K = K(T ) associated with T in the usual manner is called 
a recursive fractal. 



Figure 3. Level 2 approximations to an IFS tree and to the asso¬ 
ciated fractal. Here F = {.F(2), -F(3)} contains the IFSs generating 
SG(2) and SG( 3) respectively. Edges of the tree with a given ini¬ 
tial node are enumerated from left to right; they correspond to 
subcells enumerated anticlockwise from the bottom left corner of 
the cell corresponding to the given node. 
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Level 1 


Environment applied 


• o« ■ ■ •oon#■ •■ •o■ 



• on ■ 


Tree and prefractal approximations 


Level 0 


□ 



V=4 with types {1,2,3,4} represented by • O □ ■ 


Figure 4. Approximations to a 4-variable tree and the prefractal 
approximations to the corresponding 4-variable fractal. The IFSs 
are F( 2) and F( 3) generating SG(2) and SG( 3). The environment 
at each level is applied to the approximation at the previous level. 
The IFS labels are not shown since in this case they are determined 
by the branching number. 
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Notation 2.2. The boundary dT of a tree T is the set of infinite paths through T 
beginning at 0. 

For i £ T the cylinder set [i] C dT is the set of all infinite paths w £ dT such 
that i is an initial segment of w , written i -< w. 

The concatenation of two sequences i and j, where i is of finite length, is denoted 
by the juxtaposition ij. 

The truncation of i to the first n places is defined by i\n = i\ ... i n . 

A cut for the tree T is a finite set Ac T with the property that for every w £ dT 
there is exactly one i £ A such that i -< w. Equivalently, {[i] : i £ A} is a partition 
of dT. 


For a tree T and a node i £ T, there will usually be associated quantities such 
as an IFS a type r l £ {1,..., V} (see Definition 2.5) or a branching number 
N z . In this case i is shown as a superscript. 

In particular, the transfer operator a t acts on T to produce the tree a l T, where, 
writing for the address of node j, 

(4) :=T ij . 

That is, <7 1 T is the subtree of T which has its base (or root) node at i. 

We frequently need to multiply a sequence of quantities, or compose a sequence 
of functions, along a finite branch corresponding to a node i = i\. ..i n of T. In 
this case, i is shown as a subscript. For example, if i = i\ .. .i n then, with some 
abuse of notation for the second term, 


(5) 


li := ii, ■ ■ ■ •' := 


(F v 

l 


pF' 1 








is the product of scaling factors corresponding to the edges along the branch i\ ... i n 
and analogously for other scaling factors. Similarly, 


( 6 ) 


1 ° ° V’in : = V’zq ° VC 1 ° VC 1 * 2 ° ■ ■ ■ ° VC 1 


^12 ri 3 

is the composition of functions along the same branch. 


Notation 2.3 (Cells and Complexes). The recursive fractal K = K(T ) generated 
by T satisfies 

AT 0 

(7) K(T) = |J iff (KtfT)) = (J V’i(tfVT)), 

i= 1 \i\=n 

where the second equality comes from iterating the first. 

For |i| = n the n-complex and n-cell with address i are respectively 

(8) Ki = i^K^T)), A; := ifi(V 0 ), 

recalling that Vo is the set of essential fixed points of F £ F and is the same for 
all F. 


Assumption 2.4. In Section 3 and subsequently we assume 

Vo is the set of vertices of an equilateral tetrahedron in for some d > 2, 
Eo is the set of edges, Go = (Vo, Eq) is the complete graph on Vo- 

We will need various sequences of graph approximations {G n }^L 0 to the fractal 
K(T). In particular we use the notation G n = (V n , E n ), where 

(10) V n := (J ^i(Vo) = (J Ai, E n := (J ^(E 0 ). 

|i|=n | i | =77. |i|=n 

We can recover the fractal itself as K(T) = cliff n V n ), where cl denotes closure. 
We will write x ~ n y for x, y £ V n if x , y are connected by an edge in E n . 
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2.4. ^-Variable Trees and l 7 -Variable Fractals. Fix a natural number V. For 
motivation see Figure 4. 

The following definition of a V- variable tree and V- variable fractal is equivalent 
to that in [7] and [8], but avoids working with U-tuples of trees and fractals. 

Definition 2.5. A V-variable tree corresponding to F is an IFS construction tree 
T corresponding to F. with a type t x £ {1,...,V} associated to each node i. 
Moreover, if two nodes % and j at the same level |z| = \j\ have the same type 
r* = r J , then: 

(1) i and j have the same associated IFS F z = F J and hence the same branch¬ 
ing number N z = 

(2) comparable successor nodes ip and jp , where 1 < p < N z = N*, have the 
same type T zp = r Jp . 

The recursive fractal K = K(T) associated to a U-variablc tree T as above is 
called a V-variable fractal corresponding to F. 

The class of V-variable trees and class of V-variable fractals corresponding to 
F are denoted by fly = Uy and ICy = K.y respectively. 

Remark 2.6. A V-variable tree has at most V distinct IFSs associated to the 
nodes at each fixed level. If two nodes at the same level of a U-variable tree have 
the same type then the subtrees rooted at these two nodes are identical, i.e. 

(11) |*| = \j\ & t* = r j =» a z T = a j T. 

In particular, for each level, there are at most V distinct subtrees rooted at that 
level. 

A 1-variable tree is essentially the same as an IFS tree which generates a scale 
irregular or homogeneous fractal as in [17], [19] and [2]. 

The following is used in the construction and analysis of U-variable fractals. 

Definition 2.7. An environment E assigns to each type v £ {1,..., V} both an 

I F I 

IFS F v £ F and a sequence of types , where \F V \ is the number of functions 

in F v . We write 

(12) E=(E(1),...,E(V)), E(v) = (F v e ,t* 1 ,...,t* IF v1 ). 

For a pictorial example see Figure 4. For the following consider the case n = 2 
in Figure 4. 


Construction 2.8. A V-variable tree is constructed from a sequence of environ¬ 
ments {E k )k >i in the natural way as follows: 

Stage 0: Begin with the root node 0 and an initial type r 0 assigned to this node. 
Stage 1 : Use E 1 and the type r 0 in the natural way to assign an IFS to the level 0 
node, construct the level 1 nodes and assign a type to each of them. 

More precisely, use E 1 ^) where 


F(v) = ( 


F 


E 1 


5 ' W,1 ) • 1 • ! 


yIE® 1 1 


to assign the IFS F 0 := Fff to the node 0 and in particular determine the 
branching number V 0 := |Fff | at 0, and to assign the type t 7 := t E j to 
each level 1 node j. 


Stage n: (By the completion of stage n — 1 for n > 2, an IFS F 1 will have been 
assigned to each node i of level |i| < n — 2, all nodes i of level \j\ < n — 1 
will have been constructed and a type r 3 will have been assigned to each.) 
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Use E n in the natural way to assign an IFS to each level n — 1 node 
according to its type, to construct the level n nodes and to assign a type 
to them. 

More precisely, use E n (v) for 1 < v < V where 


E n 



E n E n 

T v,l ) ■ ■ • > T v,\Fff n | 


to assign the IFS F x := to each level n — 1 node i of type v and in 
particular to determine the branching number N l := \F^ | at the node i, 

and to assign the type r 1 - 7 := to the level n node ij. 

It follows by an easy induction that the properties in Definition 2.5 hold at all 
nodes. □ 


We now note the following facts about the connectivity properties of U-variable 
fractals. 

Lemma 2.9. Let K be a V-variable fractal. Then 

(1) K is connected 

(2) K is nested: For all i,j £ T, if [i] fl [j] = 0, then Ki fl Kj = tpi(Vo) D ipj(Vo). 

Proof. (1) The connectedness is clear as all the affine nested fractals in the family 
are connected. 

(2) In our setting this is straightforward to see as if [i] fl \j] = 0, there exists a 
k of maximal length with k -< i and k -< j, such that Ki C A4 and Kj C K 
If we write i = ki\... and j = kj i..., then i\ ^ ji and by the nesting axiom 
for F k we have K] cil fl Kkj 1 = V’fcii(K)) Fl 14j, (Vo). If the intersection is empty 
we are done. Otherwise, by our technical assumption on affine nested fractals that 
(Vo) O if? (Vo)| < 1, there is a single intersection point which is the image of 
a fixed point in Vo. If Ki fl Kj ^ 0, this is the intersection point of Ki fl Kj and 
therefore of ipi(Vo) 0 V^(Vo) as required. If Ki fl Kj = 0 we are done. □ 

2.5. Random U-Variable Trees and Random U-Variable Fractals. 


Definition 2.10. Fix a probability distribution P on F. This induces a probability 
distribution Py on the set of environments as follows. Choose the IFSs Ffj for v £ 
{1, ..., V} in an i.i.d. manner according to P. Choose types tJt for 1 < j < \F^\ in 
an i.i.d. manner according to the uniform distribution on {1,..., V} and otherwise 
independently of the F®. 

Definition 2.11. The probability distribution on the set fly of U-variable trees is 
obtained by choosing r® £ {1,...,U} according to the uniform distribution and 
independently choosing the environments at each stage in an i.i.d. manner accord¬ 
ing to Py. This probability distribution on U-variable trees induces a probability 
distribution on the set JCy of U-variable fractals. Both the probability distribution 
on trees and that on fractals are denoted by Py. We will write Ey for expectation 
with respect to Py. 

Random V-variable trees and random V-variable fractals are random labelled 
trees and random compact subsets of respectively, having the distribution Py. 
Later, when we add additional scale factors for resistance and weights associated 
with each F £ F, we will assume they are measurable with respect to F £ F. 

Although the distribution Py on environments is a product measure, this is far 
from the case for the corresponding distribution Py on fly and K,y. There is a high 
degree of dependency between the types (and hence the IFSs) assigned to different 
nodes at the same level. 
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Remark 2.12. The classes /Cy interpolate between the class of homogeneous frac¬ 
tals in the case V = 1 and the class of recursive fractals as V — > oo. The probability 
spaces (/Cy, -Py) interpolate between the natural probability distribution on homo¬ 
geneous fractals in the case V = 1 and the natural probability distribution on the 
class of recursive fractals as V —> oo. See [6] and [7]. 

Notation 2.13. It will often be convenient to identify the sample space for random 
quantities such as trees, fractals, functions associated to a branch of a tree, etc., 
with the set f ly of ^-variable trees. We use ui to denote a generic element of fly 
and combine this with other notations in the natural manner. Thus we may write 
T w , K u , iff etc. 

In particular, cr*w is the transfer operator defined in Notation 2.2 for a tree T. 
See for example the first equality in (28). However, we usually suppress ui as in the 
second equation in (28). Also see (90) and the explanation which follows it. 

2.6. Necks. The notion of a neck is critical for the analysis that follows. 

Definition 2.14. The environment E in Definition 2.7 is a neck if all are equal. 

A neck for a Y-variable tree ui is a natural number n such that the environment 
E applied at stage n in the construction of a; is a neck environment. In this case 
we say a neck occurs at level n. If i is a node in u> and |i| = n, then i is called a 
level n neck node. 

If a neck occurs at level n then the type assigned to every node at that level 
is the same. See Figure 5. It follows from Remark 2.6 that all subtrees rooted at 
level n will be the same. Note that the subtrees themselves are only constructed at 
later stages, and even the common value of the IFS at a level n neck node is not 
determined until stage n + 1. 

There is however no restriction on the IFSs occurring in a neck environment E. 
For a level n neck these IFSs are applied at level n — 1. 

Because there is an upper bound on the number of functions N F in any IFS 
F € F, there is only a finite number of type choices to be made in selecting an 
environment. It follows that necks occur infinitely often almost surely with respect 
to the probability Py defined in Definition 2.11. The sequence of neck levels in the 
construction of a Y-variablc tree or fractal is denoted by 

(13) 0 = n(0) < n( 1) < • • • < n{k) < • • • . 

The sequence of times between necks is a sequence of independent geometric 
random variables, and in particular the expected first neck satisfies 

(14) Eyn{ 1) < oo. 

Many of our future estimates rely on various a.s. properties of necks. However, 
some estimates just require that there be an infinite sequence of necks. For this 
reason we make the definition: 

(15) fly C Dy is the set of V-variable trees with an infinite sequence of necks. 

We next give an elementary result on the asymptotic behaviour of a sequence of 
geometric random variables (Ffc)fe>i- It follows that Yfc grows at most logarithmi¬ 
cally in k : and powers of Yfc grow at most geometrically, with similar results for the 
maximum and the mean of {Yi,..., Yfc}. 

The following is standard but included for completeness. Note that the Yfc need 
not actually be geometric random variables. 
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Figure 5. Compare with Figure 4, except that now a neck occurs 
at level 2. All subtrees rooted at this level will be the same, al¬ 
though they have not yet been constructed. All 2-complcxes will be 
the same up to scaling by factors determined by the construction 
up to this level. 


Lemma 2.15. Suppose {Yk}/fL 1 is a sequence of not necessarily independent ran¬ 
dom variables with P(Yk > x) < Ap x , where 0 < p < 1 and for all x > 0. Suppose 
n > 1 is a natural number. Then a.s. 


(16) 

(17) 


. Y n k 

hm sup --- < 

k —Kx) log k 


1 

logi jp 


lim sup 

k—y oo 


logfc 


< 


1 

logi !v 


lim sup 

k—>oo 


E nk -xr 

i=l zi 

k log k 


< 


2 n 

logi !v 


Proof. The case n > 1 is a direct consequence of the case n = 1, which we establish. 


Suppose e > 0. Since P(Yfc > x) < Ap x for x > 0, 


n> 


(1 + e) log k 
logi /p 


<aJ2p 


(1+e) log fc/(log 1/p) 


k >1 N ' k >1 

Hence by the first Borel-Cantelli lemma, 


A k ( 1+e l < oo. 

fc> l 


lim sup -—— < 

fc->oo log k 


1 + e 
log 1/p 


a.s. 


Since e > 0 is arbitrary, the first inequality in (16) follows. 

The second inequality in (16) is now an elementary consequence. Suppose 6 > 0. 
Using the first inequality in (16) to get the second inequality below, P a.s. there 
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exists ko = ko(u, 6) such that k > ko implies 


Yi Yi 1 

max - -- < max ;- : < 


k 0 <i<k logfc k 0 <i<klogi log 1 /p 


6. 


Hence 


Y t 1 

limsup max ;-— < -——I- <5 a.s. 


k —f oo k 0 <i<k logfc “ log 1/p 
Replacing ko by 1 and letting S —> 0 in the above implies the second inequality 
in (16). 

For (17) fix 7 > 0. Then 

Yp( fy- > lklogk \ <YYp(y > 7log E 

b \h 1 w/pj-hh io g i/J 

< ^2 kA P 7 los k/ log 1/p = A J2 kl ~ 7 < °o, if 7 > 2. 

fc>l k> 1 

By the first Borel-Cantelli lemma, if 7 > 2, 


J2i=l Y i / 7 

limsup 1 < t- 77 - a.s. 

k —Kx) fclOgfc lOg 1/p 


This gives (17). 


□ 


We also include a decomposition of sums of products of scale factors. 

It may help to note that the factors on the right side of (19) in the next Lemma 
are calculated by first choosing and fixing, for each j = 1 ... k, an arbitrary node 
of T at level n(j — 1). For fixed j all subtrees of T rooted at this level are identical 
by the definition of a neck. The factor in (19) is the sum, of products of s p type 
weights, along all paths in such a subtree starting from its root node and ending 
at a first neck level node. There is a one-one correspondence between the set of 
such paths in the subtree and the set of paths in the original tree starting from the 
chosen node at level n(j — 1 ) and ending at a level n(j) node. 


Lemma 2.16. Let Si = sf £ K for i = 1, ..., N F be scaling factors associated with 
each family F, where 


(18) 


0 < si rl f := inf{sf : i £ 1,..., N F , F £ F}, 

s sup : = supjsf : i £ 1,..., N f , F £ F} < oo. 


. (j) 

Then, writing Si = ■ ... • Si n for i = i\ .. .i n £ T, and with s\ defined in the 

natural way in the body of the proof we have 

(w) e «=n( e -?- 1 ’ 

ieT,|i|=ri(fc) 1=1 V \i\=n{j)—n(j—l) 


Moreover, 

( 20 ) 


lim — log 

k—y oo k 


E 

\i\—n(k) 


Si = Ey log ^ Si Py a.s. 


|i|=n(l) 


Proof. Let T^ denote the unique subtree of T rooted at the neck level n{k), so 
that in particular T = T. 

Then, as explained subsequently (and following the notation of (5) but with the 
F there suppressed), 


E 

i£T,\i\=n(k) 


E - 

i£T ,\i\=n(k) 


iKra(fc)- 1 ) 

M*0 


— 
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£ ..{(■ 


i^T : \i\—n(k) 


( s 0 • 


. s ili2 . 


b %2 

S i3 

i|n(l) 

i\ 
• s. 

|n(l)+l 

'*"( i)+i 


i(l) + 2 


i|(n(l)-l)' 

*"(1) j 

i|(n(2)-l) 


1 ( 2 ) 


).... 


(A 


E 

ieT(0),|i|=n(l) 

E 


x E 

iGT( fe_1 ) ,\i\—n(k)— n(k— 1) 

( 21 ) 

7 o (0) X 7 


i\n{k— 1) 

. *|n(fc- 

-1)+1 _ 

. s i|(n(fc)-lA 

Mfc—i)+i 


l)+2 

^n(fe) J 

(O),0 

S- • s 

(0),ii . £ 

,(0),iii2 

<■ ' . . 

. s (0),i|(rj(l) 

1 

*2 

^3 


J 1 )’ 0 


„(1),*1*2 

(l),i|(nl 


s i 2 

' S i 3 

• ■ • ' **„(2)_n( 

H) 







(fe-l),iii 2 


’il 

' S i2 

' S i 3 


)} 




= E 4° 

\ |i|=n(l) 


E - 

|z|=n(2)—n(l) 


( 1 ) 


E -■ 

, |i|=n(/e)— n{k— 1) 


(fc-i) 


The first and last equality are immediate from the definitions. The second equal¬ 
ity is just a bracketing of terms. 

For the third equality note that each n(j) is a neck. A term such as 
which corresponds to the edge in T from i\n(l) = i%... i n ( i) to i\ ... *n(i)*n(i)+i> is 
independent of i\n(l) and can also be regarded as corresponding to the level one 
edge from 0 to i„(i)+i of the unique tree T W rooted at every level n(l) node. Thus 
we rewrite as s^’ , with an abuse of notation in that i and i„(i)+i in the 

first term refer to words from T = T whereas i\ in the second term is the first 
element of a word from . Similarly, is also independent of i|n(l) and 

can also be regarded as corresponding to a level two edge from T^\ etc. Now use 
simple algebra to put the summations inside the parentheses. 

The final equality is a rewriting of the previous line and provides the definition 
of s) j) . 

For the Py almost sure convergence in (20) let 

= log ^ st X \ k>l. 

\i\=n(k)—n(k—l) 

By construction the Xk are i.i.d. and in particular X\ = log |*|=ri.(i) Si • By 
bounds on Sj we have 

E V \X\ \ < J2P(n(l) = n) max 11 log (A^ p s" up ) |, | log (A£ f s£ f ) |} 

n> 1 

= max 1 1 log (fV SU p 5 SUp)|5 I log (-/V in fS in f)| j-EVn(l) < 00. 

Hence, using (19), the Py almost sure convergence follows from the strong law 
of large numbers for the sequence {AA}. □ 

2.7. HausdorfF and Box Dimensions. Assume that the family of IFSs F satisfies 
the open set condition as in Section 2.2. We do not here require the affine nested 
condition. Recall the notation from Section 2.2 and from Notation 2.2. 

Splitting up and treating the necks in the manner here was done first by Scealy 
in his PhD thesis [39]. 
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Theorem 2.17. Suppose K is the random V-variable fractal generated from F. 
Then the Hausdorff and box dimension of K is Py a.s. given by the unique a such 
that 


( 22 ) 


E v log Y, V = °- 
l*l="(l) 


Proof. See the Main Theorem in Section 4.4 of [8]. The expression there for the 
pressure function is equal to the simpler expression here. This in turn leads to a 
simpler proof of that theorem, still along the lines of Lemma 5.7 in [8] but working 
with a single neck as in the (somewhat more complicated) proofs of Theorems 3.11 
and 4.14. □ 

We give a slight refinement of this result. 

Theorem 2.18. There exists a constant C such that 

1 


(23) 


lim sup - _ 

k —loo V&doglogfc 


F log h a = C, P v a.s. 

\i\=n(k) 


Proof. We can apply Lemma 2.16 with Si = li 

Since E v log Xqi|=„(i) = 0, lim^oo \ logE|i|= n (fc) = 0 a - s - Using the 

bounds (2) on N F and £ F it is easy to check that Ey ^log X)|i|=n(i) < oo. 

The law of the iterated logarithm for the sequence of random variables Xf~ = 
lo sE|i|= n{k)-n(k-i)^i now inTplies the result. □ 

3. Analysis on V- Variable Fractals 

3.1. Overview. Our V-variable affine nested gaskets are connected and nested 
by Lemma 2.9 but they need not have spatial symmetry, in contrast to the scale 
irregular nested gaskets considered in [2]. 

In order to study analysis on these V-variable affine nested fractals we define 
in Section 3.2 their Dirichlet forms and show that these are resistance forms. We 
also show that the resistance metric between points is comparable to an appropriate 
product of resistance factors. In Section 3.3 we introduce general families of weights 
and measures and prove a few basic properties. We introduce in Section 3.4 the 
notion of the cut set A&, where each cut is at a neck level and the crossing time for 
the corresponding neck cell is of order e~ k . Asymptotic properties of various quan¬ 
tities associated with these neck cells are established. In Section 3.5 we show the 
Hausdorff dimension in the resistance metric is given by the zero of an appropriate 
pressure function. 

3.2. Dirichlet and Resistance Form. The construction of the Dirichlet form 
follows [31]. 

Assume as given a harmonic structure ( D,p F ) for each IFS F in the family F. 
Since all our affine nested fractals are based on the same triangle or d-dimensional 
tetrahedron with vertices Vo, the matrix D will be independent of F and is given 

by 

(24) D(x, y) = 1, Vx, y £ Vo with x ^ y, D(x, x) = —d Vx £ Vo. 

Vectors p F = (pf ,... ,p^ F ), specifying the conductance scaling factors to be ap¬ 
plied to each cell, will be chosen to respect the symmetries of the limiting fractal. 
Assume 

1 < Pinf ■= inf(pf : 1 < i < N F , F £ F}, 

Psu P := supjpf : 1 < i < N F , F £ F} < oo. 


( 25 ) 
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The associated renormalization map for each F € F is assumed to have the usual 
fixed point property. We now state this more formally. 

Let 

( 26 ) £o{f,g) = \ E (/O)- g(y)) 

x,v&Vo 

be the Dirichlet form on the graph Go = (Vq^Eq) with conductances determined 
by the matrix D. Each edge is summed over twice, and hence the factor 1/2. 

The choice of p F is such that 
n f 


( 

(27) So(f, /) = inf < ^pf£ 0 (/io^f,/io^f) 

^ i= 1 


/i : El —> R, h\v 0 — f 


One can also regard this as placing conductors pf on each edge of the 1-cell with 
address i, which ensures that the effective resistances between vertices from Go in 
the graph Gi is the same as the effective resistance in Go itself — see Notation 2.3 
and (32). 

Define T n := {/ | / : V n —> R}. Use recursion on n > 1 to define 


AT- 


Snif,9) = E ft C(/»r,9°f) V/,ff g T\ 


(28) 


2 — 1 
iV 0 


i.e. £„(f,g) =^Pi£n-i{f °^tg o,t Pi) Vfig&Fn- 


2=1 


It follows that 
(29) 


Sn(f,g) = E Pi£o{f 9° 

|i|=n 

where pi and tpi are as in Notation 2.2. 

The sequence of forms (£ n ,F n ) can be thought of as corresponding to conduc¬ 
tances pi on the edges of the cell Aj in the graph G n , where |i| = n. 

One next defines a resistance form first on U* := (Jn>o ^ anc ^ then on its closure 
K in the standard manner as follows. By the definition of the conductance scale 
factors pf, one has monotonicity of the sequence of quadratic forms £“ (/, /). Define 

S u (f,f)= lim £“(/,/) 

n—>oo 

for f : V, —> R, restricting to those / such that the limit is finite. Using the 
definition of R in (32) with K replace by U*, one shows that R is a metric on 
U* as in Theorem 2.1.14, page 48 of [31]. Noting definition (33), one next proves 
the natural analogue of Lemma 3.2 and Corollary 3.3 for U*, without utilising 
Theorem 3.1. It follows that the metric R induces the Euclidean topology on U* 
and the completion of this metric induces the Euclidean topology on I\. 

One can now define a limit form on K by 


(30) 


£“(/,/) = lim £“(/,/) V/eJ“:= /: sup £"(/,/) < oo 

n—> oo 


where / : K —> R. Note that, from (32), if £(/, /) < oo then / is continuous and 
so is canonically determined on K by its values on the dense subset U n >o ^n C K. 

It follows from the definitions that there is a decomposition of the limit form for 
any cut A of the tree T“, see Notation 2.2. Namely, 

(31) £(f,g) = ^2 P iS a \S °ipi,g°ipi) Vf,g&F. 

i£A 
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Note the case A = {* : |i| = k} for some k and the case A = as in (52). The 
result (31) in the first of these cases with k = 1 is essentially just a consequence 
of the scaling property (28) and letting n —> oo. The general result follows from 
iterating this down the various levels corresponding to the partition A. 


The effective resistance metric between any pair of points x,y £ K is defined by 


(32) 


R(x,y) 1 = inf {<?(/,/) : f{x) = 0 ,f(y) = lj 


The proof this is a metric is essentially as in Theorem 2.1.14, page 48 of [31]. 


Recall that (£,F) is a local regular Dirichlet form on L 2 (K,y) if it has the 
following properties: 

(1) closed: F is a Hilbert space under the inner product (/, g) i—» £(f,g) + 

f fg dy, 

(2) Markov or Dirichlet: £(f,f) < £(f,f) if / is obtained by truncating / 
above by 1 and below by 0; 

(3) core or regular: if C(K) is the space of continuous functions on K then 
C(K) (~l F is dense in F in the Hilbert space sense and dense in C(K) in 
the sup norm; 

(4) local: £(f 1 g) = 0 if / and g have disjoint supports. 

For (£ , IF) to be a resistance form it is sufficient that in addition R defines a metric, 
and in particular that R(x,y) is finite and non zero if x ^ y. 

Theorem 3.1. For each uj £ fl and each finite Borel regular measure pF on 
K u with full topological support, {£ U, ,T U ) defines a local regular Dirichlet form 
on L 2 (K U , y u ). The Dirichlet form is a resistance form with resistance metric R. 

Proof. The existence of the Dirichlet form (£, F) as the limit of an increasing 
sequence of Dirichlet forms is essentially as summarised in the first paragraph of 
Section 3.4 of [31]. See [31] Appendix B3 for a discussion of Dirichlet forms. The 
proof that the Dirichlet form is a resistance form is essentially as in Section 2.3 
of [31]. ' □ 

It will be convenient here and subsequently to work with resistance scaling fac¬ 
tors which are just the inverse of the conductance scaling factors introduced in 
Section 3.2. Thus we define 

(33) r f = n= pF 1 . 

We also note that for the resistance scale factors we have 

0 < r in f := inf{rf : i £ 1,..., N F , F £ F} = p sup _1 , 

r sup := supjrf :i£ 1,..., N F , F £ F} = p in f -1 < 1. 

Next we see that the resistance metric distance between two vertices in a cell Aj 
(see (8)) is comparable to the resistance scaling factor rj for that cell. 

Lemma 3.2. There is a constant nonrandom C\ > 0 such that if x,y £ Aj and 
x ^ y then 

c-irj < R(x,y) < n. 


(35) 
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Proof. Fix x, y and i as in the statement of the lemma. 

If f{x) = 0 and f(y) = 1, then using (31), (26), monotonicity of the limit in 
(30), and (9), 

£(/,/) > Pi£o(f °i>%) = dpi, 

where d comes from the fact there are d edges in Vq containing y. 

This gives the upper bound for R in (35). 

For the lower bound, following Notation 2.2, consider a cut A of the underlying 
tree such that j = ji • ■ • j n £ A if rj is comparable to r-j. More precisely, j £ A if 

(36) r-j < n < . 

Let V = UjeA V’i(^b) be the set of vertices corresponding to cells Aj for j £ A 
(analogous to (10)). Note that i £ A and so x,y £ V. Consider the function / 
such that f(y) = 1 and f(z) = 0 for all other z £ V, and harmonically interpolate. 
Then 

(37) £(/,/)= Y Pd £ o(f °i>j) = d Y Pj<^~Pi: 

je A, y eAj je A, ye A d r ‘ nf 

using (36), taking d as in (9), and M the maximum number of regular tetrahedra 
in W 1 with disjoint interiors that can have a common vertex. 

This gives the lower bound in (35). □ 


Corollary 3.3. There is an upper bound on the diameter of the set K in the 
resistance metric, in that there exists a nonrandom constant C such that 

(38) dianiR K := sup R(x, y) < C. 

x,y{zK 

More generally, for all i £ T, 

(39) diam/jATi:= sup R(x,y)<Cri. 


Proof. First consider points x,y £ V n (see (10)) and suppose x £ Aj C V n , y £ 
A j C V n , with |i| = \j\ = n. 

Let x 0 = yo, x k £ A^*., y k £ A^ k for k = 1,..., n, with x n = x and y n = y. By 
the triangle inequality for the metric R, 

n n 

R(x,y ) < y^fi(x fc -i,x fc ) + y^i7(2/ fc -i,2/fc)- 

k= 1 fe=1 


Since x k -i,x k £ V k and all cells are triangles or tetrahedra, if a path from x k -\ 
to x k consisting of edges from V k contains two edges from the same fc-cell then it 
can be replaced by a shorter path from x k -i to x k also consisting of edges from V k . 
It follows there is a path from x k -\ to x k consisting of at most N snp edges from V k . 
Hence 

R(xk— i j Xk) A As Up r SU p, 


from (35). Hence 


R(x,y) < 


2N sup 

(1 T SU p) 


Using the density of the vertices (J n V n in K we have the result. 
The second statement follows in the same way. 


□ 


Note that the result holds for all w £ fl. 



20 


U. FREIBERG, B.M. HAMBLY, AND JOHN E. HUTCHINSON 


3.3. Weights and Measures. We next introduce a general family of measures on 
dT (see Notation 2.2) and on the corresponding fractal set K , by using a set of 
weights (wf,... ,wf F ) defined for each F £ F with wf > 0. We do not require 

E^f = !■ 


Assume 

(40) 


0 < lOjnf := infjwf : 1 < i < N F , F £ F}, 

w sup := supjwf : 1 < i < N F , F £ F} < oo. 


Following Notation 2.2 let the weight Wi of the cell Aj (corresponding complex, 
or corresponding cylinder) be the natural product of weights along the branch given 
by the node *. That is, if |*| = n, then 


(41) 


Wi := w , 


' w. 


F' 1 



Of particular interest are weights of the form wf = ( rf) a for all F £ F and 
some fixed a > 0, in which case Wi = ri a . This example is the reason we do not 
require )T) i wf = 1, since it would not be possible to achieve the normalisation 
simultaneously for all F £ F. 


The following construction is basic, and is special to the case of V- variable 
fractals. 

Definition 3.4. Let (wf,... ,wf F ) for F £ F be a set of weights as before. For 
|*| a neck let 

7/1 • 

(42) /*:=/*([<]) -. 

^|j|=|i| Wj 

The corresponding unit mass measure /r on dT is called the unit mass measure with 
weights wf. 

The pushforward measure on K under the address map 7r : dT —> K given by 
IXLi Ki\ n = {7r(*)} is also denoted by /*. 

Note that from the definition of a neck, (42) is consistent via finite additivity 
from one level of neck to the next, it extends by addition to any complex or cylinder, 
and so by standard consistency conditions it extends to a unit mass (probability) 
measure on dT. 

We note the following simple estimates for use in the rest of this subsection and 
in Lemmas 4.3, 4.4 and 5.11. 


Lemma 3.5. Suppose i and j are two nodes of the same type with |*| 
Then 


(43) 


( Win A 
\^sup / 


n 

ftj — Mi — 


/ ^sup \ 

\ Winf J 


n 


f^j • 


If i is a neck node then 


(44) 


Winf 


N f 


sup u-'sup 


n 

< IM < 



f W in f \ 

\ ^sup / 


-1 


< 1. 


|j| = n. 


Proof. Suppose N is the first neck > n. Then 


= 


E|ifc|=IV, Wik _ E|fc|=JV-n, 
E|r)!= 




|p|=A,pGT W P 


Eipi= 


|p|=Ar, P eT w p 


where wf. is the product of weights along any branch of T a of length N — n 
beginning at 0, or equivalently any branch of T of length N — n beginning at *. A 
similar expression is obtained for /ij. Since i and j are of the same type and level, 
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the trees T a and T a are identical, and so pi/wi = pj/wj. Then (43) follows from 
<„f < m < < up . 

If n is a neck then 

— 

Also we have that 


S|r>| = 


> 


W ini 


— jyn ,„n 


|p|=n W P -‘’sup‘"sup 


pi — 


]\p\—n ^ P 


< 


i + V — 

^ 'll) • 


jp|=n,p#i 


. m 


< l + (JV£f-l) 


Winf 


< 1 + 


Winf 


^sup 


W s 


n\ —1 


n\ —1 


< i, 


completing the proof of (44). 


□ 


We show in Lemma 3.7 that the pushforward measure on K is given by a similar 
expression to that for p on dT. For this we first show that the measure p on K is 
nonatomic. 


Lemma 3.6. Py a.s. we have for i £ dT 

H(i) = 0. 

Proof. Since [i\n] is a decreasing sequence of sets, from (42) 

p(i) = lim Pi\n(k)- 

k— s-oo 

By (19) and (44), writing f := iVi n f/w sup , the sequence of random variables 


(j) _ Pi\n(j) _ 

\ ^ — 
Pi\n(j-1) 2^\j\=nU)-n(j-l) W 0 




Taking logs and applying the law of large numbers we see, Py a.s., for all i, 

1 i k 

Inn -log W | n(fc) <- lim - V log (l + C"Cj)-»(j-iA = -E v log (l + C (1) ) . 

rC—^OO rC fc —>00 rv \ / \ / 


i=i 


Now, using the fact that n(l) is a geometric random variable, f < 1 and log(l+x) > 
x/2 for x < 1, we conclude 

c 


-K log (i + C n(1) ) <-\e v C {1) = 


< 0. 


2 Eyn{l)(l — C) + C 

In particular, Py almost surely, for all i £ dT , we have lim^oo Pi\ n {k) = 0 (in fact 
exponentially fast) as required. □ 


Lemma 3.7. The address map n : dT -A K is one-one except on a countable set. 
The pushforward measure p on K is nonatomic. Moreover, for \i\ a neck, 

(45) pi = p(Ki) = ———-. 

£|il=|i| u, i 

Proof. Suppose a = tt (i) = i r(j). Then for some n we have i\.. .i n = j\ ... j n 
and i n +i ^ j n + 1 - It follows that a £ K^...i n i n+1 D Ki 1 ...i n j n+1 . From Lemma 2.9, 
a £ ifii-.-inin+i (Vo)- This establishes countability of the set of points in K with more 
than one address. From Lemma 3.6 it follows this set has p -measure zero. The result 
(45) now follows from (42) and the definition of the pushforward measure. □ 
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It follows from (45) that 

(46) [ fdfj, = m [ f o fa dfi ^, 

J Ki J K° % 

where as usual p = /W is the measure on K = K u but here restricted to Ki = KV, 
and p° l = /A J,cr * is the measure on K a * = which is essentially just a scaled 

copy of the subfractal Ki. By construction, the left integral is a multiple of the 
right integral, with constant independent of /. Setting / = 1 gives the constant. 
Note that \i\ = n need not be a neck. 

The inner product (or any integral) can be decomposed as follows: 

(47) (/, g)^ = ^2Pi (/ ° ipi, 9 ° 

ie A 

for any cut A, see Notation 2.2. 

Note that (47) is analogous to the decomposition (31) for the Dirichlet form. 
The difference is that the scaling factors pi in (31) are simply computed from the 
prescribed quantities pf, unlike the scaling factors pi in (47) which are related to 
the prescribed quantities wf in a simple manner only in the case where the i are 
all neck nodes. 

We write 

(48) \\fh = (f,f) 1 J 2 

for the natural norm on L 2 (K, p). 

3.4. Time and Neck Cuts. We now introduce the special cut sets which will be 
essential for our analysis. The idea is to cut at neck nodes in such a manner that 
crossing times are comparable. 

Define 

(49) ti = pin. 

From the Einstein relation ti can be thought of as a crossing time for the continuous 
time random walk on the cell A*, with resistance given by Tj and expected jump 
time given by pi. 

Note that whereas Wi defined in (41) is a simple product of factors, as are £i, pi 
and ri following the notation of (5), this is not the case for pi and hence not for ti. 

Define 

/'cn'i Dnf Il^inf 

(50) V = N w 

and note that 0 < rj < 1. Then from (49) and (44), 

(51) p n <ti < r" up if |*| = n is a neck. 

The second inequality is clearly true for any i , not necessarily at a neck. 

Recalling from (13) the notation n(£) for the £th neck, define the cut sets of T 

(52) A o = {0}, A fe = G T : 3£^|z| = n(£), U < e~ k < j if k > 1, 

where 0 is the root node. Thus A*, is the set of neck nodes for which the crossing 
times of the corresponding cells are comparable to e~ k . 

For any i such that |z| is a neck, and in particular if % G A*,, then we define 

(53) £[i) := i if \i\ = n{£). 

That is, £(i) is the number of the neck corresponding to i. 
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We introduce further notation to capture the scale factors. 

(54) M k = |Afc|, t k = M- 1 Y, U, T k = tr 1 ; 

ieA fc 

(55) y k (i) = n(£) - n(£ - 1) if i £ A k and |i| = n(£), y k = ma xy k (i); 

i£A k 

(56) z k = max{|*| : i £ A*,}. 

Thus M k is the cardinality of the cut set A k , t k is the average crossing time for 
cells Ki with i £ or equivalently the average time scaling when passing from 
K to Ki, conversely T k is the average time scaling when passing from Ki to K for 
i £ Afc; y k (i) is the number of generations between Aj and its most recent ancestor 
also at a neck level, and y k is the maximum such number of ancestral generations 
over i £ A k ; z k is the maximum branch length of nodes in A k . 

Trivially, 

(57) min ti < t k < maxfj. 

iGAfc iGAfc 


For functions f{k ) and g(k) we will use the notation 

(58) f(k)4g(k) iff limsupF^<l. 

fc->oo 9 \k) 

That is, f{k) =4 g{k) means / is asymptotically dominated by g. 

In the next lemma we use Lemma 2.15 to estimate the asymptotic behaviour 
of y k and z k , and of the fluctuations of £(i) and ti for i £ A k . Note that sharper 
estimates for the simple case V = 1 are given in Lemma 3.9. 

Lemma 3.8. Suppose rj is as in (50). 

(a) There exist ci,C 2 > 0 such that Py a.s., if i £ A k then 

Ci/c(logfc ) -1 £{i) < C 2 k. 

(b) There exist 03,04 > 0 such that Py a.s. 

1 <yk 4 c 3 logfc, z k =4 c 4 k. 

(c) There exists /3' > 0 such that Py a.s., if i £ A k then 

k~ p 'e~ k 4 V Vk e~ k <U< e~ k . 


Proof, (a) Suppose i £ A k and let £ = £(i). From (51) and the definition of A k , 

(59) r] nW <ti< e~ k < U< r” u ( p _1) . 

In particular, n(£ — 1) < fc/ log(l/r sup ). 

It follows that 

£ = 1 + (£ ~ 1) < 1 + n(£ - 1) < 1 + --^- < c 2 k. 

log l/T S up 


On the other hand from (59), n(£) > /c/log(l/ 77 ). Using also logfc > 
log 1 / 02 , it follows from Lemma 2.15 (17), since n(£) = 1 ( n (*) — n (* — 

sum of geometric random variables, that a.s. (where i £ A k ) 

k . .. n(£)logl /?7 21 ogl/r; 


lim sup 

k—too 


< 


£{i) log k ll ^,™ P ^(log^ + log l/c 2 ) log 1/p 


=: 1 /ci. 


log^ + 

1 )) is a 


Here p is the constant probability of not obtaining a neck at any particular level > 1. 


(b) Trivially, y k > 1. By definition 

y k = max.y k (i) = max (n(£(i)) - n(£(i) - l) ) < max (n(j)-n(j- 1 )), 

ieAfc ieAfc \ / l<j<C2k \ J 
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where the inequality comes from (a). 

By Lemma 2.15 (16) with Y) = n(j) — n(j — 1), Py a.s. 

Vk , maxi cj <C2 kYj 1 

hm sup --- < hm sup ;- . -— < 


k—> oo log c 2 k - log c 2 log 1/p 


=■ c 3 , 


k— J-oo ^Og k 

where p is as in (a). 

It follows that with i £ Ak and £ = £(i), Pv a.s. 

|i| = n{£) = n{£ — 1) + n(£) — n{£ — 1) < n{£ — 1) + yk 
=4 k/ log(l/r sup ) + c 3 log k 4 k/ log(l/r sup ). 

This gives the last inequality in (b). 


(c) The third inequality in (c) is immediate from the definition of Ak- 
For the second inequality suppose i £ Ak with |i| = n(l). Then 

by a similar argument to that for the first inequality in (51). More precisely, note 
that by definition pi is a product of Aiiin(^-i) with factors that depend only on 
weights w defined along edges in the subtree rooted at i\n{£ — 1 ), followed by a 
normalisation that depends only on the same weights since |i| = n(£) is a neck. 
Hence 

U > Uw-> e~ k, q Vk 

by the definition of yk and A k- This gives the second inequality in (c). 

For the first inequality take any e > 0, in which case by (b), a.s. there exists 
k 0 = fco(w) such that k > fco implies yk < (c 3 + e) log/c, and so k > ko implies 

^Vk > ^(c 3 +e)logfc _ p-(c 3 +e) log 1/rj _ 

where /3’ = (c 3 + e) log 1/rj. Since e > 0 is arbitrary, this completes the proof. □ 
If V = 1 the above can be sharpened to the following. 

Lemma 3.9. In the case V = 1 we have the following. 

(a) There exist c\,c 2 >0 such that if i £ Ak then 

cik < £{T) < c 2 k. 

(b) There exists c 3 > 0 such that 

Vk = 1, z k < c 3 k. 

(c) There exists C 4 > 0 such that if i £ Ak then 

ae~ k <ti< e~ k . 

Proof. The first claim follows from (51) and the fact that for V = 1 every level is 
a neck. The second and third follow similarly. □ 

3.5. The Haudorff Dimension in the Resistance Metric. 

Definition 3.10. The a-dimensional Hausdorff measure of K using the resistance 
metric R is denoted by The Hausdorff dimension of K in the resistance 

metric is denoted by elf = d r ^(K). 

The following theorem is the analogue of Theorem 2.17. However, the resistance 
metric R does not scale in the same way as the standard metric in and so 
the proof needs to be modified. The proof combines ideas from Section 2 of [28], 
Section 2 of [29] and Section 4 of [ 8 ]. In the case of [ 8 ] the corresponding argu¬ 
ment is simplified here because of the use of necks. Note that we do not expect 
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the appropriate Hausdorff measure function to be a power function, unlike in [28] 
and [29]. 

Theorem 3.11. The Hausdorff dimension in the resistance metric d r j: of K is the 
unique power ao such that 

(60) E V log £ rf° = 0, 

|i|=n(l) 

Proof. This will follow from Lemmas 3.13 and 3.16. □ 

Lemma 3.12. The function 

(61) 7(a) :=E V log y rf, 

|i|=n(l) 

is finite, strictly decreasing and Lipschitz, with derivative in the interval 
[(logr inf ).EVn(l), (log T sup )E v n( 1)]. 

Since 7(0) > 0 there is a unique ao such that j(ao) = 0 and moreover ao > 0. 
Proof If a < /3, then from (34), 

7 (a) + (/3- a)(logr in f)£Vn(l) < j(fi) < 7 (a) + {/3 - a)(logr sup )E v n(l). 

This gives the Lipschitz estimate. 

Since 7(0) = Ey (log ff{i G T | |z| = n(l)}), it follows that 0 < 7(0) < 00 . 

The rest of the lemma now follows. □ 


Lemma 3.13. Suppose ao is as in Lemma 3.12. Then d r f(K) < ao, Py a.s. 
Proof Suppose a > ao- Using Corollary 3.3, 

K= \J Hi, J2 diam a R Ki<C a y rf. 

\i\=n(k) |i|=n(/c) \i\=n(k) 

From (20) and Lemma 3.12, 

lim \ log Y rf = E v log rf < 0, Py a.s. 

fc—>00 K zz ' 

\i\—n(k) |i|=n(l) 

Hence Py a.s., 

lim log rf = —00, lim rf = 0. 

k—foo ' ^ k—f 00 ' ^ 

\i\—n(k) \i\=n(k) 

Hence 'H^(itT) = 0, and so d r j(K) < ao- 


□ 


Definition 3.14. Suppose e > 0. Then A e is the cut set of T consisting of those 
nodes j = ji ■ ■ ■ j n such that 

(62) r d < e < 


Lemma 3.15. There exist non random constants c and Mi, such that for any e > 0 
and x € K, 

# {3 e A e : B ce (x) n Kj ± 0} < Ml, 


(63) 

where 


B ce (x) = {y <G K : R(x,y) < ce}. 
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Proof. Suppose x £ Ki where i € A € . 
First note 


(64) # {j eA e :KiH Kj ^ 0} < M(d + 1), 

where d+ 1 is the number of vertices of a regular tetrahedron in (recall (9)) and 
M is as in (37). This follows immediately from Lemma 2.9. 


Let V f = UjgA e V’j'(^b) =: UjeA e denote the set of vertices corresponding to 
the partition A e . 

Define u : V e —> M by u(y) = 1 if y £ Vi and u(y) = 0 otherwise. Extend u to 
u : K — > R. by harmonic extension on each Kj for j £ A £ . Then u is constant on 
Kj if Ki fl Kj = 0, and so 


£{u) 


Pj£o(u o ipj ) < M(d+ l)dpj < 

{j-.KinKjjtiH} 


Md(d + 1) 

Hnfe 


where M(d+1) is from (64) and d is the number of edges in ilj (Ao) with one vertex 
in Vi. 

Setting c = r- ln f/2Md(d + 1), it follows R(x, y) > ce if y G Kj where j £ A e and 
Ki n Kj = 0. That is, 

(65) B ce (x)nKj ^Q)=^Ki(lKj ^ 0. 

Combining (65) and (64) gives (63). □ 


Lemma 3.16. Suppose a < «o- Let p, be the unit mass measure on K constructed 
as in Definition 3.f and Lemma 3.7, with weights wf = ( rf) a for F £ F. Then 
P v a.s., for any x £ K and S > 0, p{Bs{xf) < C\8 a , where the random constant 
Ci depends on uj but not on x or 8. 

In particular, by the mass distribution principle, d r f(K) > a Py a.s., and so 
d r j(K) > «o Pv a - s - 

Proof. Fix x £ K and 8 > 0. If k is a level in the construction of T, let s(k) denote 
the first neck level > k. All balls are with respect to the resistance metric. 


From Lemma 3.15 applied to the cut Ag/ C , and with c and Mi as in that lemma, 
there are at most Mi sets Kj which meet Bg(x) and satisfy j £ Ag/ C . That is, 
satisfy, on setting j = ji ■ ■ ■ jk, 

( 66 ) rj < 8/c < r jl ,,,j k _ 1 . 

It follows that 

(67) p{B s {x)) < h{Kj), 

j G As/c ,B,5 (x)n Jf.,-^0 

and there are at most Mi terms in the sum. For each such Kj, using Lemma 3.7, 


( 68 ) 


rt*»< E = 


< 


N f 


i {k) — k 


sup 


■ rf < 


n; 


^|i|=s(fc) r i 

i(k) — k 
sup 


y|i|=s(fc) r i 


■8 a =: 6(k)8 a . 


3 — 

|i|=s(fc) i 

Here A^ sup is an upper bound for the branching number, see (2). 

We need to estimate the numerator and denominator of 9{k) in (68). For this 
we use estimates (69) and (71). 

Until we establish (71) we allow k to be an arbitrary positive integer, not neces¬ 
sarily satisfying (66). 
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Since s(k) — k is a geometric random variable, by the same argument as in 
Lemma 3.8(b), there is a constant c\ such that s(k) — k =4 ci log A; Py a.s., and so 
there is a constant 02 ( 00 ) such that 

s(k) — k < C 2 log A: Py a.s. 


for all k > 1. Hence Py a.s., for k > 1, 

(69) N$>~ k <N%f* k 

Next let 

P = E V log E rf. 

|i|=n(l) 

Then /3 > 0 since a < ao, see Lemma 3.12. It follows by (20) that as k —> oo 

i log E p v a.s. 

\i\=n(k) 

Hence for some eo = eo(w) > 0, 


(70) E rf > e 0 e fc/3/2 for k > 1, Py a.s. 

|i|=n(fc) 


However, we need an estimate similar to (70) involving s(k) rather than n(k). 
First note, by setting F) = n(i) — n(i— 1) and n = 1 in (17), that for some C 3 = 03 ( 00 ) 
we have n(k) < C 3 k\ogk Py a.s. if k > 1. Hence 


E r i^ e 0 ex P 

|i|=n(fe) 


n(k)(3 \ 
2c 3 log k ) 


for k > 1, Py a.s. 


Since n(k ) is an arbitrary neck, 


E rf > e 0 exp 
l*l=*(fc) 


s(k)p \ 
2 c 3 log k ) 


for k > 1, Py a.s. 


where k is the number of the neck s(k). Note s(k) > k. Also note that k < k. 
(Otherwise there are at least k +1 necks between levels 1 and s(k) inclusive, and so 
in particular s(k) > k. But then there are at least k necks between levels 1 and k 
inclusive, and so A is a neck. However that gives s(k) = fc, a contradiction). Hence 


(71) 


E r i - £ o exp 

|i|=s(fc) 


2c 3 log k ) 


for k > 1, Py a.s. 


It follows from (69), (71) and the definition of 0(k) in ( 68 ), that 0(k) —V 0 as 
k 00 . On the other hand, with k := \j\ we have from ( 66 ) that 


k := \j\ > log(c/ 8 )/ log(l/r min ) —>• 00 


uniformly for j € as <5 —> 0. From ( 68 ), (67) and the uniform bound M\ on 
the number of terms, there exists do = So(oo) > 0 such that 

(72) n(B$( x)) < S a for S < So Py a.s.. 

It now follows by the mass distribution principle that d r ^(K) > a Py a.s., and 
so d r f(K) > ao Py a.s. □ 
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4. Eigenvalue Counting Function 

4.1. Overview. In this section we consider random E-variable fractals constructed 
from essentially arbitrary resistances rf , from weights wf which determine a mea¬ 
sure fi, and from a probability measure P on F. See Sections 2.5, 3.2 and 3.3. With 
every realisation of such a random fractal there is an associated Dirichlet form and 
a Laplacian. The growth rate of the corresponding eigenvalue counting function is 
defined to be d s / 2, where d s is called the spectral exponent. We see in Theorem 4.14 
that Py-a.s. d s exists, is constant and is the zero of a pressure function constructed 
from the crossing times f j. The proof relies on estimates concerning the occurrences 
of necks and on a Dirichlet-Neumann bracketing argument, see Lemmas 3.8 and 4.8. 
Lemma 4.8 gives a result which holds for all realizations. (In the case wf = (rf)~ x 
some of the estimates can be sharpened, see Remark 4.9.) 

The natural metric on fractal sets constructed with resistances as here is the 
resistance metric. We saw in Theorem 3.11 that the Hausdorff dimension d} in this 
metric is given by the zero of a certain pressure function. A natural set of weights 
is wf = ( rf) df . The measure v constructed from this set of weights is called the 
flat measure with respect to the resistance metric. 

We see in Theorem 4.16 that d s {v)/ 2 = d}/{d} + 1). This establishes the analogue 
of Conjecture 4.6 in [30] for E-variable fractals. In Theorem 4.18 we show that for a 
fixed set of resistances rf , and for arbitrary weights wf and corresponding measure 
/i, the spectral exponent d s (fi) has a unique maximum when fi is the flat measure 
v. The spectral exponent in this case is called the spectral dimension associated 
with the given resistances. 

Finally, in the case of the flat measure v, we give in Theorem 4.19 an improved Py 
almost sure estimate for the counting function itself rather than its log asymptotics. 

4.2. Preliminaries. Following the notation of the previous section, we consider a 
fractal K = K u and write dK = Eo for the boundary of K. We fix a measure 
p. — pp on K and, together with the Dirichlet form £ = this allows one to 
define a Laplace operator = A" We will be interested in the spectrum of — A^ 
as this consists of positive eigenvalues. However, instead of working directly with 
—A^, we use a formulation of the Dirichlet and Neumann eigenvalue problems in 
terms of the Dirichlet form, see [31]. 

Recall the definition of F from (30). Let 

(73) r D ={/er: fix) =0,xe dK}, ££(/, /) = £“(/, /) for / € 

and let (-,-)m“ be the i nlier product on L 2 (K u ,pP). It follows as in Theorem 3.1 
that (£d,Fd) is a local regular Dirichlet form on L 2 (K \ dK,f£). Now X is a 
Dirichlet eigenvalue with eigenfunction u £ Ff, u 0, if 

(74) £“ d (u,v) = \(«,«)„- Vu £ Ff>. 

Similarly, X is a Neumann eigenvalue with eigenfunction u £ F^, u 0, if 

(75) £“(u,v) = X{u,v)^ Wv&F^. 

As usual, we will in future normally omit the dependence on to. 

By standard results [31] the Dirichlet Laplacian has a discrete spectrum 

(76) 0 < Ai < A 2 < ... where X n —> 00 as n —> 00 , 

and similarly for the Neumann Laplacian but with 0 = Ai. 

The Dirichlet and Neumann eigenvalue counting functions are defined by 

Ad(s) = max{i : Aj < s, A^ is a Dirichlet eigenvalue}, 

(77) 

Mn(s) = maxji : Aj < s, A^ is a Neumann eigenvalue}. 
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As usual, eigenvalues are counted according to their multiplicity. 

The following lemma implies the spectral exponent d s (g) in Definition 4.11 is at 
most 2 for any realization of our ’[/-variable fractals. It is used to prove the second 
estimate in Lemma 4.6. 

Lemma 4.1. With the same constant C as in Corollary 3.3, 

Md{s) < Cs , Vs > 0. 

Proof. The effective resistance between x £ K and the boundary set dK = Vq is 
defined by 

R{x,dK)~ l = inf {£(/,/) : / € T D , f(x) = l}. 

From Corollary 3.3 with the same constant C, and for any y £ dK, 

R(x,dK) < R(x,y) < C. 

The Green function for the Dirichlet problem in IF is a symmetric function g(x, y) 
which has g{x,y) < g{x,x) = R(x,dK). See, for example, Proposition 4.2 of [32]. 
In particular, 

g{x,y) < C 

independently of w. Moreover, from Theorem 4.5 of [32], 

\g{x,y) - g(x,z)\ < R(y, z). 

Hence g is continuous, and in particular uniformly Lipschitz continuous, in the 
resistance metric. 

It follows from Mercer’s theorem (for a proof of the theorem see the argument 
in [37] pages 344-345) that 

g(x,x) = (A ?) _1 &(xf 

i> 1 

and the series converges uniformly, where 4>i are the orthonormal eigenfunctions 
corresponding to the Dirichlet eigenvalues . Integrating with respect to x, 

i> 1 

for any s > 0 . □ 

4.3. Dirichlet-Neumann Bracketing. In this and the following sections, fix a 
set of weights wf as in Section 3.3 and let g be the corresponding measure. 

In order to deduce properties of the counting function for H-variable fractals we 
use the method of Dirichlet-Neumann bracketing. 

Let Afc be the sequence of cutsets (52). Using the notation of (9) and analogously 
to ( 10 ), define 

V k = IJlV’iiyo) : i 6 A fc }, 

(78) E k = |J {ifi(Eo) : i & A fc }, 

G k = {Vk,E k ). 

Thus Gk = (Vk,Ek) is the graph associated with the vertices 14 of the cells deter¬ 
mined by Afc. 
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Define ( £ k ,T k ) by 
(79) 

T k = {/ : K\V k -> M | Vi £ A k : / o ^ = /* on 

£ k (f,g) = ^2 Pi£ a \f ° 4>i, g ° A) for/, ff G J*. 

ieA fe 

The functions in J 7 ^ should be regarded as continuous functions on the disjoint 
union Uj g A fc Ar together with its natural direct sum topology. 

Define (£^,T^) by 

Fd = {/ 6 F k | Vi £ A fe /i|^ Q = 0, where fi is as in (79)}, 

£ k D (f,g) = £ k (f l9 ) for 

Thus J 7 }) is the restriction of J 7fc (and of J 7 ) to those functions which are zero on 
Vfe, and £p is the restricted energy functional. 

It is straightforward to see that 

(81) 4cJ fl cJ C J 7fc , £ k D c£ D c£ C £ fe . 

That is, £ is just the restriction to J 7 of the functional £ k and similarly for the 
other cases. 

Note that {£ k ,F k ) and (£}), J 7 })) are local regular Dirichlet forms on the spaces 
L 2 ( U i G A k Ki, n) and L 2 (K \ V k , p) respectively, with discrete spectra and bounded 
reproducing Dirichlet kernels, see [32]. 


Analogously to (75) and (74) we define the notion that A is an {£ k ,F k ), re¬ 
spectively (£p, J~p), eigenvalue with eigenfunction u. The corresponding counting 
functions are 


(82) 


A/]y (,s) = max{i : A i < s, \ is an (£ k ,T k ) eigenvalue}, 
Mq{s) = max{i : A, < s, A^ is an (£}), J 7 })) eigenvalue}. 


In order to compare the various counting functions, first note that if A = A^ 
then the decomposition (31) with £ replaced by £ k , and the decomposition (47), 
both generalise to functions f,g £ T k . The key observation now is that if A is a 
(Neumann) (£ k ,T k ) eigenvalue with eigenfunction it, then we have for all v £ T k 
that 

(83) 22 P^ V’i) = £ k {u , v) = A (it, v)/j, = A 22 t Li ( u i’i)^ ■ 

*e A k ieA k 

If we take v to be a function supported on a complex with address i £ A*,, we 
see that 

(84) £ a (u oifi,vo ifi) = ti\(u oifi,vo tpi) ,, 

since tj = pf 1 Pi- Thus ti\ is an eigenvalue of {£ a , T° ) with eigenfunction Ui = uo 
tpi. Conversely, from it* we can construct (Neumann) eigenfunctions and eigenvalues 
for (£ k ,T k ), since 

(85) 

— / \ f Ui{x) x £ K n mtifi(K) . , 

Ui{x) := < n is an eigenfunction with eigenvalue A. 

I u x kt a I I mt 'ip ‘1 ill ) 

Hence 

(86) M k {s) = 22 A/£(s) = 22 ^d{Us), 

i&A k iGA fc 

with the argument in the Dirichlet case being similar to that for the Neumann case. 
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Lemma 4.2. The following relationships hold for all s > 0 

y: A foitis) < Af D {s) < Af N (s) < Nx{tis), 

(87) ieA fc ieAfe 

Nd{s) £ jV)v(s) ^ Nd(s) + d +1. 

Proof. The proofs are a consequence of Dirichlet-Neumann bracketing and are 
straightforward extensions of those found in [31] Section 4.1 for the p.c.f. fractal 
case. The upper bound on the difference in the Neumann and Dirichlet counting 
functions is given by the number of vertices of Vo , which is d + 1 in our setting. □ 

4.4. Eigenvalue Estimates. As in the previous section, fix weights wf and the 
corresponding measure p. Let the random variable Af 1 denote the first Dirichlet 
eigenvalue. 


Lemma 4.3. If G is the upper bound on the diameter of K in the resistance metric 
given in Corollary 3.3 then for n > 2, 


( 88 ) 


C~ x < Af < 


d{d + l)p"u P 

p{K \ K b<n ) ’ 


V = 1 


A f < 


(d+l) 2 P: 


3 yj 2 

sup ^sup 


W\ 


inf 


where K b , n is the union of the d + 1 boundary n-complexes attached to the d + 1 
boundary vertices in Vq. 


Proof. Since the Dirichlet form is a resistance form we have for / € Td that 


\f( x ) - f(y) I 2 < R{x, y)£(f, /)• 


Since p is a probability measure and / € IFd, using Corollary 3.3 and the definition 
of R(x,y) in (32), it follows that 

ll/ll! < SU P \f( x )\ 2 < sup \f(x) — f(y)\ 2 < sup R(x,y)£(f,f) <C£(f,f). 

x£K u x,y£K x,yeK 


Hence by Rayleigh-Ritz, 


Af 


D 


inf 
SGJ- d 


£(fJ) 

11/112 


> c~ x . 


Next let f(x) = 0 for x £ Vo, f(x) = 1 for a; £ V n \ Vo, and harmonically 
interpolate. Then 

£(fj) = £ n (fj)<d(d + l)p^ p , 


[ f{x) 2 p(dx) > p,(K \ K b ^ n ). 
J K 


Again by Rayleigh-Ritz, this gives the upper bound. 

If V = 1 note that with n = 2 there are at least d{d + 1) interior cells as well as 
d + 1 boundary cells. Since all cells have the same type, from (43) in Lemma 3.5 
with C = w in f/w sup , 

p(K \ K ba ) > d( 2 p(K\ 2 ) = dC 2 (l - P(K \ K bt2 )) 


.-. p(K\K b , 2 ) > 


d( 2 


> 


d( 2 


1 + d( 2 ~ d + 1 ‘ 

This now gives the result for V = 1. 


□ 


In order to obtain Py-almost sure results we need to estimate the tail of the 
bottom eigenvalue random variable. Note that this result is only relevant in the 
case where V > 1 as if V = 1 then Af* is bounded above by (88). 
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Lemma 4.4. There exist constants A > 0, /3 > 0 and 7 > 0, such that 
(89) Py{ Af* > x) < Aexp(— fix 1 ). 


Proof. Let n = n u be the first level such that the following is true: if v is the 
type of a boundary n-complex of maximum mass, then at least one interior (i.e. 
non-boundary) n-complex is also of type v. 

For any such n it follows from (43) that 


w r C n C n 

if C : = — — then p(K \ K b , n ) > ——-p(A' 6;Tl ) = ——- (l - fi(K \ K b ^ n )). 


w s 


Hence 


fi(K \ K b n ) > 


d+r' ’ ' d +1 

c n 

1 + (d + l)C -n "dp 2 


> 


From Lemma 4.3 it follows that 


Af < d(d + l)(d + 2)p n s C n = d(d + l)(d + 2)f\ where £ := 


AsupU' ! si]p 

Winf 


Hence 


P v (Af > x) < P v (d(d + 1) (d + 2)C > x) 


Py n > 


log 


d{d+l)(d+2) 

log£ 


Since there are d + 1 boundary n-complexes and at least (d + l) n — (d + 1) interior 
n-complexes, and since the type of each n-complex is selected independently of each 
other n-complex, 

Py(n > y) <p( d+1 ) ! '~ 1 -( d + 1 ) =Hp (d+1) ”" 1 = A exp exp(plog(d + 1))^ 

where p is the probability that a particular complex is not of type v and A = p~ d ~ l . 
Hence setting y = log (x/(d(d + l)(d + 2)))/ log£, 

/ loffn / r \ log ( d+1 )/ log£ \ 

JV(Af > x) < A exp (^— ( ^ + 1){d + 2) ) J — Aexp(—0x y ), 

where 7 = log(d+ l)/log£ and /? = log(l/p)/(d(d + l)(d+ 2 )) 7 . □ 


Define 

(90) Ap = max{Ap ’ D : i € A*,}, 

where Xf ’ D is the first Dirichlet eigenvalue of the Dirichlet form (£^, JF) = P u ) 

with respect to the measure u . Note that A^ = Af*. 

If V = 1 by (88) we have A * < (d + l) 2 /9s Up w sup/ w inf for all k. 

Lemma 4.5. IfV >1, then with f3 and 7 as in Lemma 4 . 4 , we have Py a.s. that 

(91) A^/r 1/7 (log k) 1 /' 1 . 

Proof. In order to apply the growth estimate in the previous lemma and use 
Lemma 2.15 we use two additional properties: 

(1) The number of distinct subtrees, and hence eigenvalues, corresponding to 
each level of T is uniformly bounded (by V ); 

(2) The maximum level corresponding to nodes in Afc is asymptotically bounded 
by a multiple of k, see Lemma 3.8. 
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First consider any sequence of random ^-variable IFS trees not neces¬ 

sarily independent but all with the same distribution P = Py, see Definition 2.11. 
Let the corresponding random first eigenvalues be Yj. 

Then for all x > 0, 

P(Yj > x) < Aex p(—/3:r 7 ) = A(e~^) xl by Lemma 4.4, 

.-. P(Yy >x)< A(e~P) x , 

:. max Y? =4 /3 _1 logic Py a.s. by Lemma 2.15, 
l<j<k J 

(92) max Yj =4 /3 _1 / 7 (log k) 1 ^ Py a.s. 


For any tree T = T u there are at most V non-isomorphic subtrees rooted at each 
level. Let (Y))j>o be the sequence of random variables given by the first eigenvalue 
of T, followed by the first eigenvalues of non-isomorphic IFS subtrees of T at level 
one (there are at most V), followed by the first eigenvalues of non-isomorphic IFS 
subtrees of T at level two (again there are at most V), etc. If Yj corresponds to a 
subtree rooted at level p then by construction j < Vp. With z k as in (56) it follows 
that A^ < maxi<j<y Zfe Yj. 

Hence Py a.s., 


lim sup 

k—¥ oo 


A? 

(log k) l h 


< lim -up ln;1X| r ~' ^ ‘ 

fe-)-oo (log VZk) 7 


/ log Vz k 
\ logic 


1/7 


<r' h ( lim s«p l 2 £h±hiii) 1/1 trom (92 ) 

V fc-VOO log k J 

< r 1/T , 


since Zk =4 cjc from Lemma 3.8(b) which implies limsup i ,_ foo (log 2 fc/logic) <1. □ 


We now wish to determine the limiting behaviour of the counting function. We 
first give the following result that is true for all uj £ f V v , which we recall from (15) 
assumes that ui has an infinite number of necks. 

Recall 77 (< 1) defined in (50), and the quantities defined in (90) and (54)-(56). 

Lemma 4.6. There exists a constant Ci such that if u £ fly then 

(93) A f D (T k ) < ciM fc , M k < Af D (\*T k r,-y<) 
for all k > 0 . 

Proof. For the first estimate we have from (87), (54) and Lemma 4.1, 

A//)(Tfe) < 'y ] A/jy ( tiT k ) < (d, + 1 )M k + y ( Aff) ( t{T k ) 

iCAfe ieA fc 

< (d + 1 )M k + cT k Y,ti< c\M k . 

ieAfc 

Next note from definitions (54), (77) and (90) of M kl Ad and respectively, 
from the fact Af ' D < A 2 ' D for the equality below, and from (87) for the last 
inequality provided tf 1 < c(k) for all i £ A k , that 

(94) M k = Y, A rg (A?’ d ) < Y (Y) < Md (Af c(k)). 

ieAfc ieAfc 

But tf 1 < rj~ Vk e k < ij~ Vk T k from Lemma 3.8(c) and the definition (54) of T k . This 
gives the second estimate. □ 
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For V = 1 we can improve this to the same estimate as that obtained in [2] Sec¬ 
tion 7. 

Corollary 4.7. For V = 1 there exist constants c i and C 2 such that for all k > 0, 
M k < Afo(ciTk) and NoiTk) < c<iM k . 

Proof. As V = 1 we know that Aj 1 is bounded above and y k = 1 and hence the 
inequality on the right in (93) reduces to Mk < Nn(c\T k ) as required. □ 

We next use asymptotic information about the frequency of necks to obtain the 
following. 

Lemma 4.8. For V > 1 there exist constants C\ and a such that Py a.s. there is 
a ko(uj) for which 

A f D (T k ) < ci Af fc , M k < N D (k a T k ) if k > fcoH- 

Proof. This follows from Lemma 4.6, since Aj ^ /3 —1 /' T (log fc ) 1 / 7 by (91) and y~ yk =<: 
k& by Lemma 3.8(c). □ 

Remark 4.9. Suppose wf = pf for all F and i. Then p * = Pi/ Pj for 

i £ Afc. This measure corresponds to the random walk with constant expected 
waiting time at each node, see the sentence following (49). Note that ti = Pipf 1 = 
1/ Yjje A fc Pj is independent of i £ A fc , and that T k = J2jeA k Pj- 

For fixed k there is no spatial variability of the ti with * £ A*,, and so we can 
take c(k) = T k in (94), which improves (93) and gives 

(95) N d (T fc ) < cM k and M k <U D {\{ T fc ). 

Using (91) we can improve the second estimate to have that Py almost surely there 
is a fco(w) such that 

M k < Af D (2p- 1 ^{logk) 1 ^T k ), k > fco(w). 

4.5. Spectral Exponent. We again fix weights wf and let p be the corresponding 
measure as in Definition 3.4. 

Definition 4.10. The pressure function 7 = 7 (/?) where /3 £ R, and the constant 
fio, are defined by 

(96) 7 (/ 3 ) = £yiog ^2 7 (/ 3 o)= 0 . 

|i|=n(l) 

(It follows from Lemma 4.12 that /3q is unique.) 

The pressure function and its zero can be found computationally. See [ 8 ] for 
similar computations for the fractal dimension. 

Definition 4.11. The spectral exponent d s (p) for p is defined by 

(97) = lim 1 

2 t-t 00 log t 

We see in Theorem 4.14 that a.s. the spectral exponent exists and equals the 
constant /3 q- By Lemma 4.2 we could replace Nd by Nn- 


Recall the definition of 77 in (50) and the estimate for ti from (51). 

Lemma 4.12. The function is finite, strictly decreasing and Lipschitz, with 
derivative in the interval 


log (p 1/2 ) E v n(l),log (fsup) E v n(l) . Since 7 ( 0 ) > 0 


there is a unique j3o such that 7(/?o) = 0 and moreover /3q > 0 . 
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Proof. If a < ft then from (50) and (51), 

7(a) + ^ 2 01 (fogri)E v n{l) < 7 (/ 3 ) < 7 (a) + ^ 2 ** (log r sup )E v n(l). 

This gives the Lipschitz estimate. 

Since 7 ( 0 ) = Ey (log ff{i £ T | |i| = n(l)}), it follows that 0 < 7 ( 0 ) < 00. 

The rest of the lemma follows. □ 

Proposition 4.13. Py a.s. we have 

(98) ;lim 7 log J2 t i /2 = 7(P)- 

\i\=n(k ) 

Proof. The idea is that from the definition of a neck, log XZ |*|=n, (k) ti ' 2 is the de¬ 
ference of two random variables, each of which is the sum of k iid random variables 
having the same distribution as logX)| i |= n (i)(^« ; i ) /3/2 an d (/3/2) log X)|i|=n(i) 
respectively. 


More precisely, suppose |z| = n(k) and in particular is a neck. Then 

TiWi 


ti — 




\j\=n{k) W 3 


and so 
(99) 


( 100 ) 


i°g t i /2 = l °s ( r i“i)^ /2 -1 io § Wi ■ 

\i\=n(k) \i\=n(k) \i\—n(k) 

If we let sf = (rf wfY^ 2 or sf = wf , it follows from (25) and (40) that 
0 < Si n f := inf {sf" : i £ 1,..., N F , F £ F}, 

s sup := supjsf : z £ 1,... , N F , F £ F} < 00 . 


and we can apply Lemma 2.16. Thus (20) applied to each term on the right hand 
side of (99) gives the result. □ 

Subsequently we write Af for AI'd ■ But note that from the second line in Lemma 4.2 
the main estimates in the rest of the paper also apply immediately to ATn- 


The proof of the following theorem relies on the Dirichlet-Neumann bracketing 
result in Lemma 4.8 and the estimates in Lemma 3.8(c). 


Theorem 4.14. The spectral exponent is given by in that 


( 101 ) 


ds(fJ-) , = lim logA/~(t) 
2 t-¥ 00 logt 


§0 

2 ’ 


Py a.s. 


Proof of Theorem. Define the unit mass measure up on dT by setting, for any ft 
and for |z| = n(fc), 


v p[*\ 


J3/2 

i 

Z^\i\=n(k) 


It is straightforward to check that up is just the unit mass measure with weights 
(rjWj )^/ 2 as in Definition 3.4. 

If j(ft) < 0 or equivalently ft > fto, then from (98) for e > 0 small enough we 
have Py a.s. that there is a constant fcp such that 


up[i] > t^ 2 e > cffd 2 if fc > fc 0 . 
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As A fe is a cut set, by using the lower estimate above we have from Lemma 3.8(c) 
that Py a.s. if fc > fco then 

1 = Y oH > ct i /2 > 2 e -^/ 2 , for some /3' > 0. 

iGAfc iGAfc 

Thus 

(102) M k =$ ck pp ' /2 e kp/2 , Py a.s. 

Suppose f > 1 and let k be such that e fc_1 < t < e k . Then t < T k by 
Lemma 3.8(c) and so 

log J\f(t) logJV(T k ) log (cM fe ) /3 

—;—— < —i—7— < —j -r - =* tt, Pv a.s., 

log t log t k — 1 2 

where the second inequality is from the first estimate in Lemma 4.8 and the third 
inequality is from (102). 

As this holds for all f3 > f3 o we have 


(103) 


logA/~(s) A) 
logs ^ 2 ’ 


Pv a.s. 


Similarly we have an asymptotic lower bound. For this choose /3 < A), or 
equivalently such that 7(/3) > 0. Then for small enough e > 0 we have Py a.s. 
that for some fco = fco(w) 

vp(i) < dPJ 2 if fc > fc 0 , 
and hence from Lemma 3.8(c), that Py a.s. then 

1 = Y2 v p{P) A Yj c ^ 2 — c M k e ~ k>3 ^ 2 if fc > fco- 
iGAfc iGAfc 

Thus Py a.s. 

(104) M k > ce kf3/2 if fc > k 0 . 


From the second estimate in Lemma 4.8 and using (104) 
(105) 


log A C{k a T k ) ^ logM fc ^ P n 

fc “ fc ^ ~2 Fv& - S - 


Again choosing fc such that e k 1 < t < e k , we have from Lemma 3.8(c) that for 
some a', 

k a T k =4 k a e k < e(l + logf) a t , Py a.s. 

Hence 

\ogN(k a T k ) logAfc(e(l + logf) a 'f) 

limmf-;-- < inn ml- 1 — --, Py a.s. 

fc-> oo fc t->oo logf 

Setting y = y(t) = e(l + log t) a t, since lim^oo log y(t )/ log t = 1 and y(t) —> oo as 
t —> oo, it follows 

Iim inf < lim inf lo f W) , a.s. 

fc—^oo k *->• OO log t 

Combining this with (105), since /3 < /3q is arbitrary, implies 

log A/" (s) v A) D 

(106) —--^ TT) -?V a.s. 

logs 2 


The required result follows from (103) and (106). 


□ 
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4.6. Spectral Dimension. 

Definition 4.15. The flat measure with respect to the resistance metric is the unit 
mass measure v with weights wf = {rf) df , where d r j, is the Hausdorff dimension in 
the resistance metric (see Definition 3.4). 

The spectral dimension d s is the spectral exponent for the flat measure. 


Further justification for the definition of d s is given in Theorem 4.18. 

Recall from Theorem 3.11 that d r j is uniquely characterised by 

(107) Eylog ]T rf f = 0, 

l*l=n(l) 

As a consequence, the following theorem establishes the analogue of Conjecture 4.6 
in [30] for V- variable fractals. 

Theorem 4.16. The spectral exponent for the flat measure v is given Py a.s. by 

(' 108 ') = d / 

2 d r f + 1 

Proof. From Definition 4.15, (96), (49) and (42), if |i| = n{t) is a neck then 


nwi 


1 +d> 


S|.7l = 


d r f ' 


(109) ti := riVi = 

^-dil ="(<)' j 

Hence the spectral exponent d s {v) is uniquely characterised by 

( 110 ) 


0 = "l{ds(v)) := E v \og ^2 t ds{l/) 12 = E v log ^ 


1 +d r f 


da{y)/2 


|i|=n(l) 


|i|=n(l) \Sh-|= 


\j\=n(i) j 


Using (107), 


(i+d r f )d a {v )/2 d s {v) 


Ey log ^2 r j f = Ev lo § X! 


(l+d r f )d s {y) j 2 


0 = Ey log ^ r i 

l*l=n(l) 111=1 l*l="(l) 

Using (107) again and the uniqueness of d r j, it follows that d r j = (1 + d r j)d s {v)/ 2, 
which gives (108). □ 


We next show that the spectral dimension maximises the spectral exponent d s (p) 
over all measures p defined from a set of weights wf as in Section 3.3. A related 
result for deterministic fractals is established in Theorem A2 of [35] using Lagrange 
multipliers. Here we need a different argument, but this also establishes uniqueness 
of the wf and hence of p. 

The proof is partly motivated by [27], in particular Section 4 and the discussion 
following Corollary 2.7. We first need the following general inequality. 

Proposition 4.17. Suppose {p±,... ,pn} and {qr,..., qjy} are sets of positive real 
valued random variables, each with the same random cardinality N, on a probability 
space (f2, P). Suppose E log Pk = 0 and that the constant 7 satisfies 0 < 7 < 1. 

Then 

N / JV \ 7 

(HI) Elog^Pfc?fc < Elog ( J , 

fc=1 \fc=i / 

with equality iff q\ = • • • = qjy a.s. 
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Proof. For any TV, a suitable version of Holder’s inequality for sequences yields 

N / N \ 1-7 / N \ 7 


( 112 ) 


E Pfc «k 


k—1 


< (E^ 

\fe=i / 


\k= 1 


Taking logs and expectations, and using the assumption on the random sets 
{pi,...,Pn} gives 
(113) 

N N /NX ' 1 / N \ 7 

Ev\og^2,p k ql < (1 - 7)ElogEPfc +El°g ( ^PkQk ) = Elo S ( ) ' 

k -1 fc=1 


\k =1 


Wc=l 


This gives (111). 

If qi = •••= qw = c a.s. where c is a random variable, then equality holds in 
(111) since both sides equal Elogc 7 . 

If it is not the case that qi = ■ ■ ■ = qn a.s. then strict inequality holds in (112) 
with positive probability and hence strict inequality holds in (113). □ 

Theorem 4.18. The spectral dimension d s is the maximum spectral exponent d s (p) 
over all measures p defined from weights wf. Equality holds if and only if for some 
constant c, wf = c(rf) d f a.s., in which case the corresponding measure p is the 
flat measure with respect to the resistance metric. 

Proof. For |i| = ra(l) let pi = r ^ f , so that E v log 5 D|*|=ti(i) Pi = °- 

Suppose w = {wf | f £ F, 1 < j < N F } is a set of weights and consider the 

corresponding Wi. Let qi = Wi/r i f . 

Then from (111), 

Ev i°g E ( r i / ) 7y; i - -^Vlog^ E w ^j ■ 


|i|=n(l) 


|i|=n(l) 


Choosing 7 so that the powers of 77 and Wi are equal, gives 7 = dj/(dJ- + 1), i.e. 
7 = d s / 2. Hence 


E v log E ^ 2 = Ev lo S 


E|i|=„(i)(^^i) da/2 


|i|=n(l) 


(E|i|=n(l) W i ) 


d B /2 


< 0 . 


Moreover, by Proposition 4.17 equality holds if and only if a.s. it is the case that 
Wi/r is independent of i for |i| = n(l). Clearly, this is true iff wf = c(rf) d f a.s. 
for some constant c. 

From the definition (96) of d s (p ), we have Ey log E|ii=n(i) = 0. From 

Lemma 4.12 and the previous inequality, it follows that d s (p) < d s , and equality 
holds iff wf = c(rf) d f a.s. for some constant c □ 

We next give a sharpening of Theorem 4.14 in the case of the flat measure with 
respect to the resistance metric. This shows that for this measure, for all V > 1, 
we have the same fluctuations as observed in the version of the V = 1 case treated 
in [2], For this, let 


(114) 


<b(s) = \/s log logs, 

(j)(t) = exp ($(log t)) = exp (x/log t log log log t). 


Theorem 4.19. Suppose p is the flat measure in the resistance metric. Then there 
exist positive (non-random) constants ci, 02 , 03 , 04 , and there exists a positive finite 
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random variable cq = Co(w), such that if t > cq then 



Proof. Consider the unit mass measure 1/3 constructed in the proof of Theorem 4.14, 
where now f3 = d s is the spectral dimension as in (110). 

In the following the constant c may change from line to line, and even from one 
inequality to the next. 

If |i| is a neck and |i| = n(£) then from (108) and (109), 



(115) 


Using the law of the iterated logarithm, as in Theorem 2.18 and from the decom¬ 
position (21), Py a.s. there exists a constant c such that, for £ sufficiently large, 



Since Vd a is a unit measure and is a cut set, it follows from (115) and (116) 
by summing over i £ A*, that, for k sufficiently large, 



(117) 


ieAfc 


ieAfc 


where £{i) is defined in (53). But from Lemma 3.8(c) and Lemma 3.8(a) respec¬ 
tively, the following hold Py a.s. for i £ A k and k sufficiently large: 


c 1 k 13 e k <ti < e k , t(i) < C 2 k. 


Moreover, 4>(cfc) < c*$(fc) for some c* = c*(c) and all k > 3. It follows from (117) 
that, for k sufficiently large, 


c _ 1 M fe e _fcd<,/ 2 e _c *( fc) < 1 < cM k e~ kda/2 e c ^ k \ 


since k P d ^ 2 can be absorbed into e with a new c. That is 

(118) c -i e kd s / 2 e - c *(k) <M k < ce kd ° /2 e c * {k) . 


Given t > 0 choose k so e k 1 < t < e k . Note also from Lemma 3.8(c) that 
e fc < Tfc < cY e k , for k sufficiently large. Then from Lemma 4.8 and (118), 

(119) A f[t) < Af(T k ) < cM k < ce Ms/ 2 e c<S>(k) < ct ds/2^ t y^ 


where for the last inequality we note that 4>(fc) < $(1 + logt) < c$(logt). 
Similarly, again from Lemma 4.8 and (118), 


Af(k p "T k ) >M k > c -i e kd,/2 e -c<s>(k) > c -^t d ‘ l2 <\>(fY c . 


But Y Tfc < c(log t)P + P t < c*t for t > 2 and c* = c*(c,f3',(3"). It follows that 
Af(c*t) > c -1 ^ 3 / 2 ^^) -0 and so 


Af(t) > c-H d ^ 2 (j){t)- c 


( 120 ) 


if log log log t > 0 , hence if t > 16. 

The result follows from (119) and (120). 


□ 


Remark 4.20. By using the law of the iterated logarithm in the above we can 
show that the Weyl limit does not exist in that there is a constant c such that 
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5. On-Diagonal Heat Kernel Estimates 

5.1. Overview. The on-diagonal heat kernel is determined for resistance forms 
by the volume growth of balls. In [9] it is shown how volume estimates can be 
translated into heat kernel estimates in the case of non-uniform volume growth. 
We are in the same setting but will express the bounds in a slightly different way. 
As we have scale irregularity these will give rise to larger scale fluctuations than 
the fluctuations arising from the spatial irregularity. Note that we will establish 
bounds for the Neumann heat kernel and are in a setting where the measure is not 
volume doubling. 

In previous work, in the V = 1 setting of [2], using our notation in (54) and (114), 
the results obtained were that for all realizations there are non-random constants 
ci, c 2 such that 

ciM fe < p T ~i(x,x) < c 2 M k , Vx € K, k > 0, 

while using a sequence chosen according to Pi, there are non-random constants 
ci, C 2 , C 3 , C\ and a random variable C 5 £ ( 0 , 00 ) under Pi, such that 

cit -da/ 2 </>(l/f ) -C3 <p t (x,x) < c 2 t _ds/ 2 0(l/t) C4 , \/x € K, 0 < t < C 5 , Pi a.s. 

In the random recursive case (V = 00 ) with its natural flat measure, as considered 
in [ 22 ], the fluctuations were shown to be smaller in that there are fixed constants 
ci, C 2 , a > 0 and a random variable C 3 G ( 0 , 00 ) under P^ such that 

c 1 t~ ds/2 \\ogt\~ a <p t (x,x) < c 2 t~ ds/2 \logt\ a , VO < t < C 3 , p—a.e.xGK, a.s. 

We will show here that the on-diagonal heat kernel estimates for V variable 
fractals are determined by the local environment, see Theorems 5.5 and 5.8. In 
the case of the flat measure in the resistance metric, see Definition 4.15, we show 
in Theorem 5.13 that the global fluctuations are of the same order as the V = 1 
case for nested Sierpinski gaskets with uniform measure as described in [2]. In the 
case of a general class of measures we will see in Theorem 5.12 that /i-almost every 
x £ K does not have the same spectral exponent as the counting function (except 
when we choose the flat measure) and there will be a multifractal structure to the 
local heat kernel estimates in the same way as observed in [5], [21]. 

In order to transfer the fluctuations in the measure to the on-diagonal heat kernel 
we could apply a local Nash inequality, for example [33] or use [9]. However we use 
more bare hands arguments adapted from those of [2], [5] and [21] in order to keep 
the scale and spatial fluctuations separate. 

Note that in [9], [34] it is shown that, in the case of resistance forms with non- 
uniform volume growth and under assumptions which hold in our setting, there 
exists a heat kernel which is jointly continuous in ( t,x,y) € (0,oo) x K x K for 
every well. 

5.2. Upper Bound. We adapt the scaling argument given in [21] Appendix B to 
this setting. 

Firstly, recall from Theorem 3.1 and the definitions and discussions around (73), 
(79), (80), that (£,P), ( £d,Fd ), ( £ k ,P k ) and (Ed^d) are local regular Dirichlet 
forms on L 2 (K, p), L 2 (K' \U 0 , p), T 2 (U ie A fc I<i, p), L 2 (K \ V k , p) respectively. For 
A > 0 let 

£a(/, 3)=£(/,2) + A(/, 5 ) m , 

with similar expressions for the other Dirichlet forms. The space T equipped with 
norm E 1 ^ 2 is again a reproducing kernel Hilbert space and we write g\,g\, gx,g\' D 
for the corresponding reproducing kernels. 

We state a scaling property of the Dirichlet form. 
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Lemma 5.1. For all f,g £ T we have 

£\(f,9) = Pi £ xuU°^9 °4>i)- 

ieAfc 


Proof. This follows by the scaling in (31) and (47) and the definiton of ti in (49). □ 

Let g^' a be the reproducing kernel associated with the Dirichlet form £ a D A on 
K a with Dirichlet boundary conditions and let gf be the reproducing kernel for 
the Dirichlet form £f on K a with Neumann boundary conditions. 

Lemma 5.2. We have for all i £ and x £ Ki, that 

9\’ a = Pi9x/^( x , x )- 

and 

9\ (V^O). ^ i l (■ x )) = Pi9\lu ( x > x ). 


Proof We consider g^’ a ftp. 1 (x),ip i X (x)), for x £ Ki, which is the reproducing 

kernel for (<?£),a>-^d) on L 2 {K^, ). We note that g®' a {4’i 1 {y),'4>i 1 {x)) = 0 

for all y £ K\Ki. Using this, the reproducing kernel property and the scaling, we 
have for x £ Ki, 


9x;:\4f\ x ),^\ x )) 


£dM’ d M,9&\k\WZ 1 {x))) 

je A fc 

Pi£D,xtA9x D (M^ x ),9xi a \,^\x))) 

Pi9x D {x,x) 


as required. 

The second equation follows by the same argument. 


□ 


It is straightforward to see that, as 

and gx{x,x) = [inf {£x(f,f) '■ f 6 F,f{x) > l}] -1 , (with similar expressions for 
9\,9x D ,9x) we have 

(121) gx’ D (x, x) < 9x (x, X) < gx(x, x) < gx{x, x), Mx £ K\V k . 

Lemma 5.3. There exists a function C{ A) such that for all A < oo 

sup gx(x,x) < C( A) < oo. 

xeK 

Proof We follow the proof of [1] Theorem 7.20. Note that for any fixed x £ K we 
have gx(x, .) £ T and hence using (32) 

\g\(x,y) - gx(x,x)\ 2 < R(x,y)£ x (gx{x, .),g\(x ,.)). 

By the reproducing kernel property and the global bound on the resistance across 
K from Corollary 3.3 we have 

\gx(x, y) - gx{x, x)\ 2 < Cg x {x, x). 

Rearranging 

9\(x,y) > g\{x,x) - ( Cgx(x,x)) 1/2 , 
and integrating over y against /i we have 

gx(x,x) < j + [Cgx{x,x)) 1/2 . 
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The result then follows easily. □ 

Lemma 5.4. There exists a constant C such that for all i £ and x £ int(Ki), 

g t -i(x,x) < Cpf 1 . 

Proof. By Lemma 5.2 and (121) we have for x £ Ki 

9\’° = P%9x^{x,x) < pig\/ u {x,x) 

< Pi9\/ti ( x > x ) = 9\ 

Now set A = 1 and note that by Lemma 5.3 g\ is uniformly bounded. Thus 

9?' a < pig 1/t .{x,x) < gf (i>i 1 {x),ipT k {x)) < C. 

Rearranging we have 

g t -i(x,x) < Cpf 1 , 

as required. □ 


Theorem 5.5. There exists a constant c such that 

p ti (x,x) < cpf 1 , Vx £ Ki,Vi £ A fc . 

Proof. As 

/» OO 

g\(x,x)= / e~ xt p t (x,x)dt , 

Jo 

we have, by the monotonicity of pt(x, x) in t, that for all u 


f u 1 

9x(x,x) > Pu(x,x) / e~ xt dt = p u (x,x) — 
Jo 


A u 


x 

Thus, setting A = tf 1 = 1/u, we have 

p ti (x,x){ 1 - e~ 1 )ti < g t -i(x,x) < Cpf 1 . 
Rearranging and the definition of ti then gives the result. 


□ 


5.3. Lower Bound. We follow a standard approach see for instance [2], [5]. For 
this we require an estimate on the exit time distribution for balls. We start with 
some preliminary results. 

Let {X t : t > 0} be the diffusion with law P associated with the Dirichlet 
form (£, P). We write P^ for the law of the process with Xq = x and E x for the 
corresponding expectation. We write Ta = inf {t > 0 : X t G A} for the first hitting 
time of the set A. For * £ we write 

Di= { J [Kj : Kj n Ki ± 0} 

ieAfc 

for the union of the complex Ki and its neighbours. Let A k (z) := {j G A* : z G Kj}. 
For z £ 14 we define 

D\,= (J K j: dD%~ (J MV 0 )\{z}, dDi := )J 1>j(V 0 )\MVo)- 

j£A k (z) jeA k (z) 

We will also use the notation dKi := ipifVo). 

Recalling (2), (40) and (55) we let yfP = Yntk° Vu V = N infW inl /N sup w sup and 
write x(/e, no ) = {'q/fi) Vk f) Vk . 

Lemma 5.6. There exist constants Ci and no such that 

cix{k,no)e~ k < E x Tg Di < sup E z Tg Di < C 2 e~ k , \/x G Ki,Mi G A k . 

z&Ki 
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Proof. We begin by observing that 

(122) E x Tg D . = E x Tg K . + J2 FX ( X T aKi = v)WT BDi . 

yGdKi 

To treat the first term we note that the Dirichlet form restricted to Ki with 
Dirichlet boundary conditions is a reproducing kernel Hilbert space with the associ¬ 
ated Green function gx^x ,.) as the kernel. Let f(y) = gx i (x,y)/gx i (x,x). By the 
definition of / and the reproducing kernel property we have <?(/, /) = l/gxi{x,x). 
By the definition of the effective resistance we also have that gx t (x, x ) = R(x, dKi). 
As gxi is harmonic away from x and is 0 on dKi we have that 0 < f(y) < 1 for 
all y. Hence, putting these observations together and using Corollary 3.3, we have 
that, for any y £ Ki, 

E v T dK .= ( gKi(y, z )p{dz) < R{y,dKi)y,(Ki) < criHi < c e ~ k , 

J Ki 


as i £ Afc. 

We next consider the exit time from Di started at a point y £ dKi. 

Let Uq = 0 and set Ui = inf 1 1 > Ui-i : X t £ Vfc\{A'[/ i _ 1 }|. Then Xi = Xjj i is a 

discrete time Markov chain on 14- Let S = inf{n : X n £ dDi}. By construction we 
see that {X n : n < S} can be viewed as a Vo + 1 state discrete time Markov chain 
with Vo states as the vertices of Ki and an absorbing state given by amalgamating 
the vertices in dDi. By construction this Markov chain has transition probabilities 
given by the conductances on Gk- As two of the vertices in dKi must be internal 
to a triangle or d-dimensional tretrahedron in Ki\\i\—x the conductance between 
the edges across and at least one edge to dDi are comparable or otherwise the 
conductances across Aj are smaller and hence ES < oo independent of k. 

The time taken for the original process to exit is then E y Us■ We now compute 
the time for a step. 

The same argument as before for the first term in (122) but using g D y gives 

^ v T dD y = f g D y(y,z)y.(dz) < R(y,dD?)/j,(D?). 

’ Jd v 


Now observe that by the definition of resistance we have 


R{y,dD\)<R{y,z), Vz£dD y . 

Thus we have R(y,dD y ) < min z ^gx>v ^?(2 A~)- By our estimate on the resistance 
in Lemma 5.6 this gives R(y, dD y ) < minj g A fc (y) r j- Hence, as the number of cells 
that meet at y is bounded, 


E : v T bd v < 


min r.j 

j£^k{y) 


Y. kj 

je^k(y) 


< c max rjtij < ce 
' j&A k (y) 


We are now ready to show E v Tgr>. < Ce k . To see this we use 

S 

E y T dDi = E y U s = E y ^ (Ui - U^). 

i=l 

Note that S’ is a stopping time with respect to {Eu % })L 0 , where is the 

filtration generated by X. As E(Ui — L4—i |-7 r c/ i _ 1 ) = E Xui ~ 1 Ty k \{ Xu I, a minor 
modification of Wald’s identity shows that 

E y T dDi < ce~ k E v S, My £ dKi. 


Putting this back into (122) gives the upper bound. 
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For the mean hitting time lower bound we return to (122) to see that 

E x Tgjj. > nrin E v Tgr>. > min E v Tq D v. 

' y£dKi 1 yedKi <■ 

Using the properties of <?£>?, and setting f(z) = go? (y, z)/gD? {Vi 2 /), we see that 


\Hy)-f(z)\ 2 <R(y,z)£(fJ) = 


R(y,z) 


R(y,z) 


g D y ( y, y) R{y , dD \)' 


Let 


A c y := {z:R(y,z)<cR(y,dD v i )}. 


Let j* £ Ak{y) denote the index at which min 3 - eAfc ( :(/ ) rj is attained. Thus, by the 
boundedness of |Afc(x/)|, we have 

(123) R(y,dDy) > arj*. 

We now show that A y must have measure comparable with . 

By decomposing the cell Kj* we have 

i<r = U 4’rd(K(v rj T)), 

j : \j\ =n 

and we write k with |fc| = n such that y £ t/jj*k{K(&'' k T )) = Kj*k- Then, by 
Corollary 3.3, for any 2 £ Kj* k we have a constant c such that 

R(y,z) < crj* k < crj*r” up , 

and hence by (123) 

CT n 

R{y,z)<-^R{y,dD\). 

^1 

Thus, if we take no = inf{n : r™ up < ci/c} and set C 2 = — r -^ R , we have Kj* k C A y 2 
where C 2 < 1 . 

Hence for 2 £ Kj* k we have | f(y) — f(z) | 2 < C 2 - As f(y) = 1 we see that we 
must have f(z) > d = 1 — y/cf. Thus for any y £ dK\ we have, writing k no + |j*| 
for the first neck after no + \j*\, 


E V T; 


dD y 


gDv(y,z)n{dz) > c'g D y(y,y)n(K^ k ) 


ID? 


= c'R(y, dDDnj 


J2i w ki ■ \ki\ = k no ,i £ T a3 } 


3 T,{ w i : 1*1 = fcn 0 >* e T<Tr } 

Now apply (123), the upper and lower bounds on the weights and Lemma 3.8, 

Wnt 


E v T dDi > c 3 v r g r N^ no 


A^sup^sup 


> QQ-k'qVkfjyk+i? — \-yk+n 0 

as required. 

Lemma 5.7. There exist constants c 3 , C 4 such that for x £ Ki, i £ A k 
P x (T dDi <t)< 1 - c 3 x(k,n 0 ), for t < cj^x(k, n 0 ) 2 e~ k . 


□ 


Proof. We note that 


TdDi <t + I{T aD .>t}(Td Di ~ t). 


Taking expectations 


E x T aDi < t + E* (l {TaDi>t} E x *T dDi ) 

< t + ¥ x (Tgu i > t) sup E v Tgjg i . 

y£Di 
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Rearranging and then applying our exit time estimates from Lemma 5.6 

P x (T dDi <t) < - t - w — + l - EXT 7y T 

sup ygD . E y T dD . sup yg£ j. E y T aD . 

< Cl e k tx(k, no) -1 + 1 - c 2 x(fc, n 0 ). 

Thus, if t < ^c 2 c]" 1 x(fc, no) 2 e _fe , we have 

P x (T aD . <t)< 1 - ic 2 x(fc,n 0 ), 

as required. □ 

Theorem 5.8. There are constants c,a' such that for t < C 4 e~ k x(k,no) 2 

P 2 t{x,x) > cx(k, no) 2 fj.(Di)- 1 , Mx G A'j,i G A fc . 

Proof. A standard argument gives the following. If t < ^C 4 x(&, rio) 2 e~ k , then by 
Lemma 5.7 

(c 2 x(k,n 0 )) 2 < P x (X t G A) 2 = ( [ p t {x, y)p(dy)) 2 < p(Di)p 2t {x,x), 

J Di 

as required. □ 

Finally we can use the estimates on yi~ to provide a Py a.s. estimate in terms of 
the scale factors. 

Theorem 5.9. There are constants c, (3 such that Py a.s. for sufficiently large k, 
for t < ce~ k k~ 2 P 

Pt{x,x ) > cp(Di)~ l k~ 2p ,Mx G Ki,i G A k . 

Remark 5.10. In a different setting [5] obtained a finer estimate on the exit time 
from a complex which enables the derivation of a finer form of this on-diagonal 
estimate. We do not derive such a result here though we expect that the same 
techniques could be applied to do so. Our result is enough to enable us to compute 
the /i-almost everywhere local spectral exponent 

5.4. Local Spectral Exponent. As in [5] we will see that the local spectral di¬ 
mension obtained by considering the limit as k —> oo of pt t ( x , x) for x £ Ki, i G A k 
will in general not coincide with the global spectral dimension. 

We have the following preliminary result. Let i x G dT be such that iFiqfc —> {x} 
as k oo. 

Lemma 5.11. There exists a constant c such that -Diq n (fc+[ c logfc]) C A-V \ n (k) f or 
all sufficiently large k for p-a.e. x G K, Py a.s. 

Proof. Let T n ( m \ h denote the addresses of the three boundary cells at the m-th 
neck. By Lemma 3.5 we must have 

a := Ey max u, < 1. 
jeT na)ib 

Now for i G {j G T : \j\ = n(k + m)} we have Di C A' i |„( fc ) if A'j fl dKi\ n ^ k ) = 0. 
Then, setting A = {Rj : i G {j G T : |j| = n{k + to)}, Ki fl dK^ n ^ f 0}, we have 

= Ey ^ Mi^{RT i n9Rr i | n(i , ) 7 ‘0} 

ie{jeT:\j\=n(k+m)} 

= Ey J2 W E ^ 

ie{j£T:\j\=n(lc)} j£T n (k+m), b 


Eyp(A ) 
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By construction the terms 22 = are independent and equal in distribution 

to allowing us to write 

m 

EyfJ,(A) = Ey ^2 PiEy e 

ie{jeT:|j|=n(fe)} i£Tn(m),b 3 =1 

< (d + l)a m . 

Thus we have 

oo oo 

Ey ^ h(-T £ E . D^x | n (fc_{_[ c l 0 g fc]) ^ — ^1 ^ a ^ ^ 

fe=l /c=l 

for large enough c. Hence Py a.s. we have 


£ Ff . jn(fc+[c log fc]) -P-I :c |n(/c) UO.) 0, 

as required. 


□ 


For the rest of this section we write T’n(i) = {j £ T : \j\ = n(l)} for the tree 
up to the first neck. Take another set of weights {{ref }|=i}f£F satisfying the 
conditions of Section 3.3 and define the associated measure A- 
Observe that by (44) and the definition of 77 we have 

Mi>(—r (1) , ii>r? n(1) , ier„d). 

Hnf 

Thus log pi > n(l) log ^ and as 

0 > V' Ailog/^i > n(l)log—, 

Hnf 

we have 

Py| ^2 fiilogpi\ < cE v n{ 1) < 00 . 
ieT n(1) 

We can control Ey\ J] igT (1) Ail°g^*l hi the same way. 

In the same way as [5] we can now determine the local spectral exponent for the 
heat kernel pt(x,x) defined with respect to the reference measure /i for fi almost 
every x. 


Theorem 5.12. Py almost surely, for fi-almost every x £ K we have 

lim logPt^;^) _ d s (fi) _ Ey Sier n(1) A* l°g Pi 
t^O log t 2 Ey Ei G T„ (1) Ailogfi ' 

Proof. For x £ K we have a sequence i\n(k) for which Di\ n ( k ) —> {x} as k —»■ 00 . 
By monotonicity of the diagonal heat kernel in time for t £ (ii|n(fc)’ ti\ n (k— 1 )] we 
have p t (x, x) < Pt i]nW (x, x) and thus 

logp t (x,x) \ogp t (x,x) 

iimsup--- < limsup —--, Py a.s. 

t—>0 k-lg t k—t oo log ti| n (fe-l) 

Now using Theorem 5.5 we have 

k 

l °SPt i{nm (x,x) < C-logpi\ n(k) =C-Elog • 

j-l Pi\n(j-X) 

We now consider the probability measure dfidPy on {1,..., N sup } x fly (with the 
product tr-algebra). If the point x is chosen according to A, then the terms (id 1 = 
„ l|n(j) are independent and equal in distribution to Ui\ n (i) under dpdPy. We can 

also express — logfj| n (fc) in terms of independent random variables t^ 1 ' 1 defined in 
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the same way. It is easy to see that log< : i| rJ ,(fe)/ log iiintfe—i) —5► 1 for any x £ K, 
Py-almost surely and hence 


(124) 


Um l°gp,(x.x) 
t->o - log t 


< limsup 

k—¥ OO 


lQ g^|n( fc )(^^) 
log ti\ n {k) 


< limsup 

k—too 




As the mean of log fi^' ) is finite we can apply the strong law of large numbers under 
dfidPy to see that 


lim — V logM J = Ey Mi log Mi, u a.e. x £ K, Py a.s. 

k-t oo K 

i=i ier n (i) 


Similarly we can find the limit for the denominator in (124). Thus we have 
l °gPt iln(k) {x,x) ^vEieT„ (1) AilogMi 

/c—>-oo log^i|n(fc) Ey Z-/igT n (i) l^S 

For the lower bound we define £(i, k) = £ if i\n(£) £ A*,. Thus 


fog ti\n(£(i,k)—l) ^ k ^ log ti\n(£(i,k)) • 


Hence, it is clear that, by the independence 
(125) 

,. £(i,k) -£ -1 

lim —-— = lim --= ———- 

k ^>oo k I-+00 LOg ti\ n (£) wL»GT ( i) 


Ai log ti 


, A a.e. a; G AT, Py a.s. 


Now, from Theorem 5.9, we have that Py a.s. for ce f fc+1 l(fc + l) @ < t < ce k k & 
we have for x £ Ki , i £ A*,, 


logP 2 t(x,x) 
-log t 


> 


log(cp(Dj) 1 k 2/3 ) 
log (ce _ ( fc+1 l (fc + 1) _2/3 ) 
log c - 2/3 log k - log p(Dj) 
log c — k — 1 — 2/3 log (k + 1) 


Thus 


log Pt(x,o;) 
t A*6 — log t 


lim 


= lim 

t —^0 


> lim 

k—t oo 


log p 2 t{x,x) 
-log t 
-log p(Di) 


We now observe that by Lemma 5.11 we have a constant c' such that 

log i\n(J(i.k ))) — log Mi|n(K(i,fc) —[c' log K(i,fc)]) • 

Using this, (125) and writing £(i,k ) = £(i,k) — [d \og£(i, fc)], we have 


lim 

k—too 


log Mi 




= — lim 

k—too 


£(i, k ) 


fc 

lim 

k—too 


£(i,k) 


E vJ2 i( : 


i=1 

Ai log Mi 


lo g/4' 


0 ) 


E v E*eT n(1) Ai log U 


for A a.e. x £ K 1 Py almost surely, as required. 


□ 
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In the case where the reference measure p is the flat measure v in the resistance 


metric, the weights are proportional to rP and FVlogEie 
calculation shows that 


rP = 0, a simple 


d> s 

~2 


E vT,i 


i£T n (i) 


Hi log Hi 


E vJ2i 


ieT "(i) 


E vJ2i 


* eT ™( i) 


Hi log ti 
Hi log ■ 


21 j£ 


r»(l) 


E vJ2i 


ier r 


(i) 


Hi log ■ 


d f E v Eis 


Hi log n 


(1 + d r f )E v £i 6 T Hi log n 


l s 


d r f +1 2 

Indeed in this case we can go further and give a bound on the size of the scale 
fluctuations. 

Theorem 5.13. If v is the flat measure in the resistance metric we have constants 
Ci, C 2 , C 3 , C 4 £ (0, oo) and a random variable 0 < C 5 such that Py a.s. for any x £ K 

c 1 (f(l/f)~ c ' 2 t~ ds/2 < p t { x,x) < Cz(j){l/t) c H~ ds/2 , 0 <t <c 5 . 

Proof. We begin by observing that for i £ A& we have ti < e~ k and thus substitut¬ 
ing in the upper bound estimate from Theorem 5.5, for x £ Ki 

(126) p e -k(x,x) < p ti {x,x) < cvf 1 . 

d r 

By (116) we have that Py almost surely for sufficiently large k. >rP exp(—c$(^(*))) 
and hence r\ +df exp(—c$(^(i))) < ti < e~ k . Thus, using Lemma 3.8, 

n < e- fe /( 1+rf /) exp(c'i>(^(i))), 

and 

d r f 

p e -k(x,x) < e 1+d f exp(c"<!>(£(*))). 

Thus, for e~ k <t< e~ k+1 , and as l < ck < —clogt, we have for any x £ K, 

p t (x,x) < Ctr ds/2 exp(c'<F(log (1/t))) = Ct~ de ^ 2 <j)(l/t) c , Py a.s. 

For the lower bound we observe from Theorem 5.9 that Py almost surely for 
sufficiently large k, for t < ce~ k k~ 2d 

Pt(x,x) > c/i(-Di) _1 fc _2/3 , Va; £ Ki,i £ A*. 

As e~ k > tj = fjVj > r^ +df exp(—c<f>(£(j))) Py almost surely, as before we have 

vj < e- fc /(!+ rf 5) e xp(c'$(£(j))). 

Then as the number of cells in Di is bounded and £(j) < ck by Lemma 3.8(a), we 
have 

r d f 

"( D i) = - --- dfhKinx^} 


J-6A fc T,j>:\j>]= n (t(j)) V 

< ^2 r d, I {KinKj jiii,}exp(c^(£{j))) 

je A fc 

< ce- fe ^ /(1+ ^ ) exp(c"$(^(j)))- 
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Thus, Py a.s. for sufficiently large k for t < ce~ k k ~ 2 ^, 

Pt(x,x) > ck~ 2 Pe kd f^ 1+d f' > exp (— c"$(/c)), Vx € iv*. 

For ce~^ k+1 \k + 1 ) _2/3 < t < ce~ k k~ 2/3 we have C\ee~ k k~ 2/3 < t so that e k k 2 ^ > 
C 2^ 1 and 

p t (x, x) > c fc- 2 ( 2 ^+ 1 )/3/(< i ?+ 1 )i- rf J/( 1 +< i ?) exp(-c"$(fc)), Vx e A'j. 

Now as k < log c + log (1/i) we have for sufficiently small 0 < t, for any x £ K 
Pt{x,x) > b'\ logt |”' 3 t~ d ^ 2 exp(— c"<&(\ logi|). 

By adjusting c" we can absorb the logarithm into the exponential term and we have 
the result. □ 
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